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1.  INTRODUCTION 

It  is  now  generally  accepted  that  the  behaviour  of  a  reactive  material 
needs  to  be  considered  in  terms  of  its  explosive,  its  thermal  and  its  mech¬ 
anical  properties.  Most  explosives  deform  before  initiation  and  the  details 
of  this  deformation  determine  how  energy  can  be  localised  to  give  "hot  spots". 
Once  heat  is  produced  locally  the  "hot  spot"  temperature  is  controlled  by 
the  heat  being  produced  and  the  flow  of  heat  away  from  the  "hot  spot".  To 
understand  this  process  it  is  important  to  have  the  reaction  parameters  of 
the  explosive  and  thermal  properties  such  as  specific  heat  and  conductivity. 
For  many  explosives,  particularly  new  compositions,  these  are  not  always 
known  accurately. 

The  research  on  this  contract  has  considered  the  mechanical  and  thermal 
properties  in  some  detail.  Studies  of  initiation  and  propagation  have  been 
greatly  aided  by  high-speed  photography.  Research  on  reaction  kinetics  and 
thermal  properties  has  used  mass  spectroscopy,  differential  scanning  calori¬ 
metry  (DSC),  thermogravimetric  analysis  (TG)  and  scanning  electron  microscopy. 
Mechanical  properties  experiments  include  measurements  of  hardness,  yield 
strength,  coefficient  of  friction  and  fracture  surface  energy. 

The  four  areas  of  work  described  in  this  report  have  cither  been 
published  (sections  2, 3, A)  or  are  ready  for  submission  (section  I). 

The  texts  of  these  papers  are  given  in  full. 

2.  THE  IGNITION  OF  A  THIN  LAYER  OF  EXPLOSIVE  BY  IMPACT;  THE  EFFECT  OF  ADDED 
POLYMER  PARTICLES.  G.M.  Swallowc  and  J.E.  Field. 

2.1  Introduction 


A  study  of  the  deformation  behaviour  of  a  thin  layer  of  material  when 
impacted  is  relevant  to  a  range  of  problems  including  the  sensitiveness  of 
explosives.  A  standard  test  procedure  for  assessing  the  hazard  involved  in 
handling  explosive  materials  is  to  impact  a  sample  (typically  20-50  mg)  with 
a  falling  weight.  The  sensitiveness  of  a  sample  is  usually  expressed  in 
terms  of  a  50%  drop  height  to  cause  ignition  (i.e.  the  height  from  which  the 
weight  would  ignite  50%  of  a  series  of  samples).  There  are  standard  methods 
for  obtaining  the  50%  height  (see,  for  example,  Dixon  and  Massey,  1 957 ).  However 
although  the  test  has  been  used  for  many  decades  it  is  only  recently  that 
a  systematic  attempt  has  been  made  to  photograph  the  sample  during  its 
deformation  (Heavens,  1973,  Heavens  and  Field  197A).  A  combination  of 
high-speed  photography  combined  with  pressure  measuring  techniques  helped 
establish  the  physical  processes  occurring  during  ignition  and  propagation. 

In  general  it  is  thought  that  the  initiation  of  an  explosive  by  mech¬ 
anical  shock  is  thermal  in  origin,  although  some  workers  believe  that  a 
tribochemical  or  a  molecular  fracture  mechanism  may  be  responsible  in  certain 
circumstances  (Taylor  and  Weale,  1932,  UbbeLohde,  19A8,  Fox,  1970).  Oil  the 
basis  of  localized  thermal,  energy  or  "hot  spots"  as  the  source  of  the  ex¬ 
plosion,  four  possible  mechanisms  have  been  envisaged  for  ignition  by  impact. 

(i)  Adiabatic  compression  of  trapped  gas  spaces  (Bowden,  Mulcahy,  Vines 
and  Yoffe,  I9A7,.  Chaudhri  and  Field,  1974). 

(ii)  Viscous  heating  of  material  rapidly  extruded  between  the  impacting 
surfaces  (Eirich  and  Tabor  1947,  Bolkhovitinov  and  Pokhil  *958)  or  by  capi¬ 
llary  flow  between  grains  (Rideal  and  Robertson,  1948). 
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(iii)  Friction  between  the  impacting  surfaces  and/or  grit  particles 
and/or  grains  of  the  mnterinl  (Bowden  and  Gurton,  1941). 

(iv)  Localised  adiabatic  deformation  of  the  thin  layer  upon  mechanical 
failure (Afiums'ov  and  Bobo lev,  197 1 ,  Winter  and  Field,  1975). 

It  was  shown  quite  convincingly  (Bowden  et  al,  1947)  that  the  presence 
of  gas  bubblos  affected  the  impact  sensitiveness  of  nitroglycerine, and 
Bowden  maintained  repeatedly  (Bowden  and  Yoffe,  1949,  1952,  1958,  Bowden, 

1950,  1963)  that  mechanism  (i)  was  the  primary  cause  of  ignition,  not  only 
in  liquids  but  also  in  solid  explosives.  This  view  has  with  equal  persis¬ 
tence  been  rejected  by  the  Russian  school  of  thought  (Andreev,  Mnurina  and 
Ruankova  1955,  Bolkhovitinov,  1959,  Ainnns'ev  and  Bobolev, 1971 ) .  In  impact 
experiments  on  solid  explosive  materials  it  was  discovered  (Kholevo,  1946) 
that  the  sensitiveness  was  considerably  reduced  if  the  sample  was  prevented 
from  flowing.  Afanas'ev  and  Bobt>lev(l96l)  produced  evidence  from  strain- 
gauge  experiments  which  suggested  that  during  impact  the  sample  suddenly 
undergoes  a  type  of  fracture.  These  authors  proposed  a  mechanism  of  hot¬ 
spot  formation  by  the  release  of  energy  along  slip  surfaces  formed  at  the. 
onset  of  mechanical  failure  (Afanas'ev  and  Boboley, 1971,  Afanas'ev,  Bobolev, 
Kasarova  and  Karabanov  1972). 

The  high-speed  camera  work  of  Heavens  and  Field  (1974)  showed  that 
impacted  'samples  in  the  drop-weight  test  may  undergo  bulk  plastic  flow, 
show  evidence  of  partial  fusion  and  even  (with  PETN)  melt  completely.  The 
flow  speed  during  these  processes  is  considerable  and  may  reach  a  few  100  ms. 
Ignition  occurred  at  a  small  number  of  local  hot-spots,  Strain  gauge  measure¬ 
ments  showed  that  pressures  of  typically  0.5  -  1  £»Pa  (5  -  10  kbar)  were  ach¬ 
ieved,  If  o  sample  failed  by  plastic  flow  this  was  accompanied  by  a  sharp 
pressure  dmp,  It  was  nf^er  t’.ho  pressure  drop  that  rhe  high  velocity  flow 
phenomena  took  place  and  Ignition  occurred.  Hot-spots  were  thought  to  be 
caused  either  by  local  obstructions  in  the  flow  or  by  the  rapid  closure  of 
gas. pockets  trapped  in  the  material. 


Recent  work  has  been  concerned  with  studying  the  effect  of  added  grit 
particles  in  greater  detail.  The  case  of  the  hard,  high  melting  point 
particle  is  reasonably  well  understood  (Bowden  ond  Gurton  1949’  Bowden  and 
Yoffe  1952,  Ubbclholdo  1948).  The  basic  idea  is  that  vbon  two  solids  rub  or 
impact  together  the  hot-spot  temperature  at  the  interfaco  is  determined  by 
the  solid  with  the  iowost  melting  point.  Properties  such  as  thermal  con¬ 
ductivity  and  hardnesa  are  important  but  only  as  second-order  effects.  In 
fact  by  chcsing  grits  of  different  melting  point  and  measuring  the  impact 
sensitiveness  Bowden  and  Gurton  were  able  to  measure  hot-spot  ignition 
temperatures  for  a  range  of  explosives.  The  hot-spol.  temperature  found 
by  these  workers  for  TETN  in  a  drop-weight  impact  situation  was  cs.  700  K 
for  hot-spots  of  micron  dimensions. 

It  is  now  known  that  certain  polymers  such  as  polycarbonate  (PC)  and 
polysulphone  (PS)  can  also  sensitise  explosives  In  the  drop-weight  test 
though  others  such  as  polypropylene  (I'P)  do  .not  (Buun  1973).  At  first 
sight  this  result  appears  surprising  since  polymers  arc  rtla*-4 ve.ly  soft 
And  have  softening  temperatures  usually  well  below  the  700  K  mentioned 
above.  The  primary  object  of  our  new  work  wan  to  try  to  understand  rhe 
mechanisms  by  which  come  polymers  give  increased  sensitiveness  and  why  others 
do  not.  At  the  start  it  appeared  possible  that  the  explanation  could  bo  a 
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chemical  one  based  on  the  production  of  free  radical#  which  aided  the 
reaction  kinetics  of  the  explosive  or  a  mechanical  one  in  which  hot-spots 
were  produced  by  the  deformation  of  the  polymer.  The  first  experiments 
were  concerned  with  possible  chemical  effects  and  used  techniques  such  as 
differential  scanning  calorimetry  (DSC)  and  thermogrevimetric  analysis  (TO. 
The  mechanical  effects  were  studied  using  the  photographic  and  pressure 
measuring  techniques  described  earlier. 

2.2  Studies  of  Chemical  Effects 

Experimental 

These  experiments  involved  measuring  the  reaction  parameters  for  decom¬ 
position  of  pure  explosive  and  explosive/polymer  mixtures  and  assessing 
any  differences.  The  apparatus  used  was  a  Stanton-Redcrof t  TG-750  thermo¬ 
balance  which  allows  the  weight  of  the  sample  to  be  monitored  whilst  it  is 
heated  at  a  constant  rate,  thus  giving  a  weight-loss  versus  temperature 
curve.  The  samples  used  consisted  of''' I  mg  lots  of  a  mixture  containing 
25X  (by  weight)  of  polymer  to  75X  explosive.  This  percentage  of  polymer 
ie  much  greater  than  likely  to  be  found  in  practice,  but  was  choaen  to 
amplify  any  interaction.  The  explosive  used  in  all  these  experiments  was 
pentaerythritol  tatranitrata  (PETN).  The  analysis  follows  the  method 
of  Hauser  end  Field  (1978  and  Bection  3  of  this  report)  and  so  is  only 
briefly  outlined  here. 

Decomposition  is  assumed  to  follow  an  equation  of  the  form 


^  -  A  exp(-E/RT)wn  0) 

with  w  the  fractional  raaidual  weight  at  time  t,  E  the  activation  enargy, 

A  tha  frequency  factor,  n  tho  reaction  order  and  T  the  temperature.  Since 
the  heating  rate  H  is  constant,  one  may  write  dT/dt  -  H.  Substituting 
in  equation  1  and  taking  logsyieldst 

In^r  “  In"  -  |T  +  n  In  w  (2) 

Thus  a  plot  of  (In  dv/dT  -  n  In  w)  against  1/T  gives  a  straight  line  if 
n  haa  been  choaen  correctly.  The  value  of  E  is  determined  from  the  slope 
of  the  line  and  A  from  its  intercept.  The  analysis  was  carried  out  on  the 
Cambridge  IBM  370/165  computer.  The  programs  used  fit  u  cubic  epline  to 
the  original  data  and  produce  differentials  which  are  then  used  to  give  plots 
of  (ln4*  -  n  In  w)  against  l/T  for  values  of  n  from  0  to  2  in  staps  of  0.2. 
Tha  straightast  section  of  these  curves  for  each  stage  of  tha  reaction 
there  may  be  more  than  one)  was  then  selected  by  eye  end  the  appropriate 
values  of  n,  E  and  A  obtained.  The  method  is  illustrated  by  the  activation 
energy  plot  for  a  sample  of  pentaerythritol  tetranitrate  (PETN)  in  Figure  I. 
The  technique  was  assessed  in  two  ways.  Firstly  by  using  the  obtained 
reaction  parameters  and  simulating  a  decomposition  curve  which  could  be 
competed  with  the  original.  Secondly,  the  dehydration  of  calcium  oxalate 
monohydrate  was  observed. 

CaC2041120  -*•  CaC204  ♦  HjO 
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This  reaction  is  well  understood  and  the  reaction  parameters  are  well  docu¬ 
mented  in  the  literature  (Gurrieri,  Siracusa,  Cali,  1974,  Freeman,  Carroll, 
1958).  the  results  obtained  with  the  technique  outlined  above  were  in  good 
agreement  with  the  literature  values  and  are  set  out  in  table  I. 

Table  1:  Dehydration  of  Calcium  Oxalate  Monohydrate 

E/kJ  mol"1  In  A/ In  e~l 

This  work  89.5  *  3  20.2 

Literature  88.6  21.1 


Results 

(a)  Thermal  decomposition  of  PETN/polymer  staples 

The  polymers  reported  to  sensitise  PETR  (Bean,  1973)  were  polycarbon¬ 
ate  (PC),  polysulphone  amber  (PS)  and  to  a  lesser  extent  polypropylene  (PP)< 
Results  obtained  in  this  work  confirmed  the  dramatic  affect  produced  by  PS 
and  PC  but  showed  very  little  increased  sensitivity  in  the  presence  of  PP *. 
Samples  of  these  polymers  supplied  by  AWRE  were  ground  in  a  fresser  mill  to 
sizes  of  <76  wn  and  then  mixed  with  PETN  to  produce  experimental  samples. 

The  heating  rate  used  throughout  was  10  deg  min  ‘  and  the  purger  gee 
nitrogen  or  argon.  The  results  are  summarised  in  Table  2,  the  figures 
listed  being  averages  of  between  three  and  seven  different  experiments. 


Table  2;  Effect  of  Polymeric  Additives  on  the  Thermal 
~~  Decomposition oi  PETN'  " 

Sample  E§/kJ  mol*-1  Ed/kJ  mol"1  In  As/ln  s~‘  In  Ad/ln  s”1  n  T/K 


PRTN 

143 

+ 

5 

193  -  7 

3) 

+ 

4 

44 

-  6 

0.7 

478 

PETN+PC 

130 

+ 

mm 

10  * 

214  -30 

30 

+ 

6 

53 

i.o 

1.0 

476 

PETN+PP 

122 

+ 

,5,* 

172  -15 

26 

+ 

4 

41 

it 

0,6 

477 

PETN+PS 

no 

♦ 

10 

231  -20 

30 

4 
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In  each  csse^a  two  stage  reaction  was  observed  (Figure  1),  the  first 
stage  being  identified  with  the  sublimation  of  PETN  (E  ,A  )  and  the  second 
its  decomposition  (E.,A,).  The  values  obtained  for  the  activation  energies 
of  sublimation  043  *a5  i?J  mol”')and  decomposition  (193  *  7  kJ  mol"’*)  are  in 
good  agreement  vtth  the  literature  values.  Cundall,  Palmer  and  Wood, (1977) 
quote  150.4  kJ  mol~'  for  the  sublimation  of  PETN  while  Roberston  (1948), 
Rogers  end  Harris  (1966)  Maycock  end  Vernsker(1970)  and  Ng,  Field  and  Hauser 
(1976)  give  values  for  the  decomposition  of  197,198,188  and  192  kJ  mol"*' 
respectively.  The  value  for  In  Aa  of  43.5  is  also  in  good  agreement  with 
the  45.6  reported  by  Robertson  (1948). 

A  larger  spread  ia  observed  in  the  values  of  in  the  PETN/additive 
samples  than  the  pure  PETN.  However,  although  some  activation  energy 
differences  are  observed,  there  ere  no  significant  changes  or  consistent 
trends.  The  values  of  frequency  factor  are  much  less  accurate  than  those 
of  activation  energy  since  they  have  been  obtained  by  measuring  the  intercept 


*  The  sensitising  effect  reported  by  Bean  for  PP  in  PETN  is  surprising 
and  we  query  the  effect.  It  is  contrary  to  our  results  and  in  fact  Bean's 
results  for  PP  in  other  explosives  show  little  or  no  positive  effect... 


rather  than  the  slope.  Again fwirhin  the  errors  of  measurement, they  show  no 
significant  changes  from  the  values  for  pure  PETN.  If  the  additive  had 
significantly  affected  the  reaction  towards  greatly  increased  reaction  rate 
then  A  should  have  increased  and  E  decreased.  Table  2  shows  that  only  small 
changes  in  A  and  l£  took  place,  and  further  that  any  changes  of  A  and  E 
tended  to  cancel  (i.e.  as  A  increases,  E  also  increases  and  vice  versa). 

The  only  consistently  observable  feature  in  the  results  is  a  change 
in  the  value  of  n  from  0,7  to  1.0  (except  in  the  case  of  PP) .  A  value  of 
n  ca.  0,7  implies  a  surface  reaction  (for  which  n  ■  2/3)  and  n  *  1.0 
is  the  value  one  would  expect  for  a  bulk  reaction.  The  explanation  for 
this  is  that  the  high  proportion  (25X)  of  filler  material  distributed 
throughout  the  sample  breaks  up  the  PETK  giving  it  a  greater  surface  area, 
thus  favouring  a  bulk  rather  than  a  purely  surface  reaction.  The  result 
for  PP  in  explained  by  the  fact  that  it  softens  and  melts  at  a  lower  temp¬ 
erature  than  the  other  polymers.  PETN/PP  samples  were  observed  Co ‘give 
droplets  of  polymer  which  then  coalesced  to  form  one  drop.  This  effectively 
removed  the  dispersion  of  particles  which  gives  n  approaching  unity  leaving 
a  surface  reaction  with  n  ca.  2/3. 

(b)  Thermal  decomposition  of  PETN/Benzoyl  peroxide  samples 

The  results  of  the  precoeding  section  show  that  the  additives  tested 
have  very  little  effect  on  the  thermal  decomposition  of  PETN.  It  may, 
therefore  be  concluded  that  their  sensitising  effect  is  not  due  to  polymer 
decom^oiition  products  interacting  with  PETN.  However,  the  possibility 
still  exists  that,  when  polymers  fracture  or  flow  during  impact,  large  con¬ 
centrations  of  free  radicals  may  be  produced,  and  these  could  react  with 
PETN  and  cause  ignition.  Walker  and  Green  (1976)  have  proposed  that  the 
mechanical  production  of  free  radicals  is  important  in  initiation  and  have 
done  some  experiments  in  which  they  find  that  tetramethylammonium  tribo- 
hydride  (a  free  radical  donor)  causes  increased  reaction  with  ammonium 
nitrate. 


In  order  to  test  this  hypothesis,  mixtures  of  benzoyl  peroxide  and 
PETN  were  used  in  a  series  of  TG  experiments.  Benzoyl  peroxide  is  an 
organic  oxide  which  thermally  breaks  down  in  the  temperature  range  310 
to  350  K  to  form  two  identical  free  radicals 


C 
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C  -  vs 


2  C  .Hr 
6  5 


2CftH5 


+  2C02 


It  is  widely  used  in  the  polymer  industry  as  a  ’starter'  material  for 
free  radical  polymerisation  as  it  provides  a  good  source  of  easily  formed 
free  radicals, 

A  number  of  TG  runs  wore  curried  out  on  both  benroyl  peroxide  by 
itself  and  two  PETN/benzoyl  peroxide  mixtures.  The  mixtures  used  had  com¬ 
positions  of  75%  to  25%  and  93%  to  7%  by  weight  of  PETN/benzoyl  peroxide. 


Due  to  the  decomposition  of  the  benzoyl  peroxide,  the  only  activation 
energy  thnt  can  be  estimated  is  the  decomposition  of  PETN,  the  sublimation 
stage  being  buried  in  the  bcnuoyl  peroxide  reaction.  Simulated  computer 
decompositions  of  these  mixtures  were  mode  by  combining  75%  of  the  numerical 
value  of  pure  PETN  run  with  25%  of  the  numeric.il  value  of  a  pure  benzoyl 
peroxide  run  (and  similarly  for  the  93/7  mix).  The  figures  thus  obtained 


Jv. 


V 


were  treate.  if  they  had  been  real  data  and  the  analysis  carried  out  in 
the  usual  manner.  Comparison  of  the  E,A  and  n  values  obtained  in  real 
and  simulated  experiments,  Table  3,  reveals  that  what  appears  at  first  to 
be  a  very  strong  chemical  interaction  is  due  to  the  combined  effects  of 
th>.'  two  materials  acting  independently. 


Sample 


Table  3i  Reaction  Parameters  for  PETN/Benzoyl  Peroxide  (b.p.)  Samples 
E./kJ  mol  *  In  A./ln  S  '  n 


7%  b.p.  93%  PETN 


185  -  8 


46  -  2 


7%  b.p.  93%  PETN 
(Simulation) 

25%  b.p.  75%  PETN 


176  +16 


120  -  6 


4i  is 


29  -  2 


25%  b.p.  75%  PETN 
(Simulation) 


143  +  20 


33-2 


(c)  Effect  of  ultra-violet  light  on  the  thermal  decomposition  of  PETN/grit 


samples 


As  a  further  test  of  the  free-radical  initiation  hypothesis,  PETN/ 
additive  samples  were  irradiated  with  ultra-violet  light  aB  they  were  being 
heated  in  the  TG.  U-V  is  frequently  effective  in  degrading  polymers  by 
breaking  chains  with  the  formation  of  free  radicals.  Both  PP  and  PS  are 
susceptible  to  photodegradation  by  U-V  whereas  PC  is  degraded  relatively 
much  Less.  It  would  therefore  be  expected  that,  if  free  radicals  are 
important  in  sensitising  PETN,  both  PP  and  PS  in  combination  with  U-V  should 
show  an  appreciable  effect  and  PC  little  or  no  effect.  The  results  obtained 
in  these  experiments  are  set  out  in  Table  4. 

Table  4:  Reaction  Parameters  in  the  Presence  of  U-V  Light 


Samnle  E  /kJ  mol  '  In  A  /In  s  '  E„/kJ  mol  '  In  A./ln  S  *  n  T_/K 


PETN  +  U-V 


143  -  5 


120-4 


31-4 


28-2 


193  -  7 


44-6 


183  -  16  45  -  4.5 


0.7  478 


0.6  458 


PETN  +  U-V  + 
PC 


124  -  8 


29-3 


196  -25  49-7 


0.9  458 


PETN  +  U-V  + 
PP 


102  -  8 


22-4 


186  -  8 


45  ±  3 


0.6  458 


PETN  +  U-V  + 

PS 


112  -  8 


25  -  2.5  Ed, 368  -  9  InAd,  99  -  6 
Ed j 1 85  -\5  lnAd2  *7  ’  5 


1.4  446 


Except  for  the  rather  anomalous  behaviour  of  PS  in  the  decomposition  stage 
(see  later),  the  results  indicate  that  only  small  effects  are  produced. 


8f. 


The  changes  in  the  sublimation  parameters  E  and  A  tend  to  cancel  in  all 
cases  and  for  the  decomposition  parameters  and  changes  again  effect¬ 
ively  compensate,  particularly  for  polycarbonate.  PP,  PS  and  pure  PETN 

ail  show  decreases  in  Ej  (E^  in  the  case  of  PS)  combined  with  increases 
in  Ad  (Ad2  with  PS).  Thus  the  reaction  is  helped,  albeit  a  little,  by 
the  U-V  radiation,  This  is  illustrated  by  the  slight  but  consistent, 
reduction  in  the  reaction  peak  temperatures  Tp. 

PS  shows  the  most  significant  changes  in  that  an  extra  reaction  stage 
appears,  This  result  was  found  to  be  very  repeatable,  though  it  is  not 
I'u] ly  understood.  Although  Ed|  (for  the  new  stage)  is  very  high,  the  value 
of  Ad |  more  than  compensates  for  this,  and  Tp  is  reduced  by  12  K.  The 
initial  and  final  stages  of  the  reaction  are  virtually  unchanged  by  the 
presence  of  this  new  stage,  which  only  has  n  small  (ca.  9  K)  range  of  activity. 

A  further  series  of  experiments  on  the  FETN/PS  system  was  carried  out 
in  which  the  U-V  light  was  switched  on  or  off  during  the  course  of  a  run 
and  the  effect  on  the  reaction  parameters  noted.  These  experiments  show 
that  the  stage  of  reaction  at  which  the  U-V  is  present  is  quite  important. 

If  the  U-V  is  switched  on  at  any  stage  before  a  third  of  the  material  present 
has  decomposed,  then  a  three  stage  reaction  is  observed,  otherwise  no  effect 
is  observed.  When  U-V  has  been  present  for  the  initial  stages  of  the  reaction 
(even  if  only  one  eighth  of  the  explosive  has  decomposed)  extinguishing  the 
light  has  no  effect  and  the  reaction  proceeds  in  three  stages,  as  if  the 
U-V  were  present  throughout. 

These  results  indicate  that  the  U-V  is  effective  in  breaking  down  the 
PS  to  produce  free  radicals  which  then  influence  the  decomposition  of  PETN. 

If,  however,  PETN  decomposition  has  proceeded  to  a  significant  extent  before 
the  U-V  is  present,  the  irradiation  has  little  effect.  This  is  probably 
because  the  free  radical  concentration  produced  by  the  PETN  alone  is  already 
much  greater  than  that  which  the  lamp  cm  produce.  The  observation  that 
extinguishing  the  lamp  early  in  the  reaction  does  not  stop  the  three  stage 
process  indicates  that  the  major  influence  is  the  concentration  of  free 
radicals  in  the  initial  stages  of  decomposition.  The  influence  of  U-V  may 
thus  be  summarised  as  producing  such  a  concentration  of  radicals  in  the 
early  stages  of  PETN  decomposition  that  the  decomposition  process  is  speeded 
up  and  the  reaction  parameters  (and  presumably ’the  details  of  the  mechanism) 
are  changed. 

That  PC  is  not  influenced  in  this  way  is  not  surprising  since,  in 
comparison  with  other  polymers,  polycarbonates  are  not  degraded  to  any 
Large  extent  by  U-V  (Ranby  and  Rabak,  1975).  It  is,  however,  surprising 
that  if  such  an  effect  should  be  produced  with  PS,PF  remains  ineffecitve. 

Both  PS  and  PP  are  degraded  to  a  large  extent  by  U-V  and  both  produce 
large  concentrations  of  radicals  when  degraded.  It  suggests  that  the 
sulphune  groups  have  a  specific  effect  on  at  least  one  of  the  stages  of  the  PETN 
decomposition.  Support  for  this  comes  from  the  work  of  Reich  (1973)  who  tested 
the  compatibility  of  a  range  of  polymers  with  various  explosives.  PS  is  one 
of  the  polymers  he  lists  as  incompatible  with  PETN  since  he  detected  an 
increase  in  the  value  of  the  reaction  parameter  n  in  the  present  of  PS. 
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From  the  series  of  experiments  described  in  this  section,  it  may  be 
concluded  that  there  is  no  major  chemical  reaction  between  PETN  and  any 
of  the  polymer  fillers.  The  only  exception  to  this  is  PS  and  then  only  in 

the  extreme  conditions  of  having  a  high  intensity  of  U-V  shone  on  Che  sample, 
However,  the  effect,  even  with  PS,  is  still  only  marginal  and  the  conclusion 
has  to  be  that  chemical  effects,  although  present,  are  only  of  secondary 
importance  compared  with  the  mechanical  ones  to  be  discussed  in  the  following 
section. 


2 .  1  Photographic  Observations  of  the  Impact  Behaviour  of  Polymers 
and  Explosives 

(a)  High-speed  photography  of  Impacts  on  polymer  samples 

The  impact  behaviour  of  both  polymers  and  explosives  was  photographed 
using  a  C4  rotating  mirror  framing  camera  at  framing  intervals  of 

ca.  6.5  ps.  Impacts  were  observed  using  a  transparent  toughened  glass 
anvil  system  similar  to  that  of  Heavens  and  Field  (1974)  and  is  shown 
schematically  in  Figure  2.  The  weight  of  5.5  kg  was  dropped  by  an  elec¬ 
tromagnet  from  a  height  of  1  m,  the  fall  being  guided  by  three  cylindrical 
rods. 


The  behaviour  of  small  discs  of  polymer  ca.  2  ram  diameter  and  0.8  mm 
height  under  impact  is  illustrated  in  Figure  3.  Figure  3  shows  selected 
frames  from  an  impact  on  a  disc  of  PP.  A  trace  of  radial  expansion  versus 
time  is  given  in  Figure  5.  With  this  material,  there  are  no  rapid  changes 
in  the  expansion  rate  and  bulk  deformation  occurs  throughout.  Figure  4  is 
a  higher  magnification  sequence  for  an  initially  2  mm  diameter  disc  of  PS. 
Its  radial  expansion  is  also  included  in  Figure  5.  This  material  deforms 
initially  in  a  bulk  manner  but  eventually  fails  catastrophically  (see 
frame  4(c)  and  Figure  5).  PC  gives  similar  catastrophic  failure  to  PS. 
Microscopic  examination  of  the  samples  after  Impact  confirms  that  catas¬ 
trophic  failure  of  PC  and  PS  is  associated  with  rapid  cracking  and  shearing 
in  localised  bands.  Figure  6  is  a  micrograph  of  a  deformed  PS  disc.  There 
is  a  network  of  fracture  throughout  the  sample  and  sets  of  fine  parallel 
shear  bands  (one  set  is  labelled  s).  An  enlarged  view  of  region  4  ,  but 
rotated  so  that  the  bands  are  horizontal  in  the  figure  is  given  in  6b.  It 
is  argued  later  that  it  is  these  regions  of  localised  deformation  (crack¬ 
ing  and  shear  banding)  which  give  rise  to  "hot-spots". 

These  observations  lead  to  the  conclusion  that  the  polymers  used  may 
be  classified  in  two  broad  groups!  (i)  Polymers  such  as  PC  and  PS  which  fail 
by  the  production  of  many  fast  moving  cracks  and  local  shear  bands.  This 
is  accompanied  by  rapid  radial  expansion,  (ii)  Polymers  such  as  PP  which 
deform  plastically  at  quite  a  high  rate  but  do  so  by  bulk  deformation  and 
without  cracking  and/or  shear  banding. 

(b)  High-speed  photography  of  Impacts  on  explosives 


Heavens  and  Field  (1974)  photographed  the  impact  of  a  wide  range 
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of  granular  explosives.  Powdered  PETN  was  observed  to  be  compacted  into 
a  pellet  and  then  to  undergo  severe  plastic  deformation.  The  material 
became  translucent  and  finally  completely  transparent  with  very  rapdid flow 
(300  m  s-l)  which  Heavens  and  Field  attributed  to  melting.  Ignition  occurred 
in  this  final  stage,  RDX  partially  fused  wnen  impacted  between  glass  anvils 
and  there  was  indirect  evidence  of  complete  sample  melting  of  RDX  between 
steel  anvils.  Pressure-time  curves  obtained  by  Heavens  and  Field  showed 
that  a  sharp  pressure  drop  occurred  when  the  sample  failed  plastically. 
Ignition,  in  samples  without  added  grit  particles,  always  occurred  after 
the  pressure  drop  and  when  the  explosive  was  flowing  rapidly.  In  the 
present  work  a  large  number  of  experiments  were  carried  out  on  pellets  of 
PETN  and  the  behaviour  was  very  similar  to  that  described  by  Heavens 
and  Field  for  granular  material.  Figure  7  illustrates  a  typical 
result  in  which  a  25  mm  pellet  of  PETN,  pressed  to  5k  bar  and  of  diameter 
5  mm  is  impacted  by  a  5.5  kg  mass  from  a  height  of  i  m.  The  same  type  of 
behaviour  as  described  for  powders,  can  clearly  be  observed.  In  frames 
a  to  f  the  material  is  deformed  plastically  and  becomes  translucent.  Trans¬ 
parency  and  rapid  flow  (230  m  s  1)  occur  in  frame  g  with  ignition  at  sites  A 
and  B  almost  simultaneously. 

The  results  obtained  with  a  large  number  of  such  sequences  show  that 
ignition  only  occurs  if  tapid  flow,  as  observed  in  frame  g,  takes  place. 
Although  rapid  flow  is  a  necessary  prerequisite  for  ignition  the  converse 
is  not  necessarily  true  and  many  failures  occur  without  a  resulting  ignition 
This  suggests  that  there  has  to  be  a  particular  region  in  the  flowing  material 
where  a  "hot-spot"  develops.  These  regions  could  be  associated  with  gas 
bubbles,  foreign  particles  or  enhanced  flow.  If  discontinuities  of  sufficient 
size  do  not  exist,  ignition  does  not  take  place.  Note  that  fast  reaction 
always  starts  at  localised  regions  (see  Figure  7  and  the  sequences  in 
Heavens  and  Field). 

« 

(c)  High-speed  photography  of  impacts  on  explosive/polymer  samples 

A  number  of  high-speed  photographic  sequences  at  framing  rates  of  up 
to  2  x  l()5  frames  per  second  have  been  taken  of  impacts  on  explosive/grit 
samples.  These  sequences  were  taken  in  transmitted  light  using  the  experi¬ 
mental  aet-up  described  above.  The  first  point  to  emerge  from  these 
experiments  was  that, unlike  impacts  on  pure  PETN  where  ignition  never 
occurs  before  transparencv  and  rapid  flow  of  the  sample, have  occurred, 
when  a  grit  is  present,  ignition  can  take  place  before  the  explosive  has 
been  sufficiently  compacted  to  fail  plastical iy T  form  a  transparent  layer 
etc.  This  observation  implies  that  the  pressures  needed  to  cause  explosion 
are  much  smaller  than  with  a  pure  material. 

If  a  sample  ignites  while  it  is  still  opaque  it  means  that  small 
polymer  particles  within  the  pellet  can  not  be  seen,  and  thus  association 
of  the  initiation  sites  and  grit  particles  is  difficult.  To  overcome  this 
problem,  two  different  approaches  were  adopted.  In  the  first  instance, 
pressed  pellets  of  pure  explosive  were  surrounded  by  polymer  discs  and  the 
impact  phenomena  of  the  whole  array  photographed.  The  second  type  of  experi¬ 
ment  involved  pressing  pellets  of  explosive,  into  which  had  been  put  single 
polymer  discs,  effectively  modelling  a  single  large  polymer  particle  in  a 
sample  of  pure  explosive.  Examples  of  results  obtained  from  these  experi¬ 
ments  are  illustrated  in  Figures  8-11. 

Figure  8  shows  a  0,65  mm  high  PETN  pellet  surrounded  by  five  I  mm 
high  PS  discs  .  During  the  impact  the  FETN  disc  (dark  central  area)  breaks 
up  and  is  forced  into  close  contact  with  the  polymer.  The  polymer  discs 
deform  plastically  and  then  undergo  catastrophic  failure  as  described 
earlier.  In  frame  8b  all  the  discs  are  exhibiting  cracking  and  shear  band- 


ing  and  by  frame  8c  have  failed  catatrophically .  In  frame  8d  we  can 
see  the  first  ignition  site  (labelled  I)  with  reaction  spreading  along 
two  paths  into  the  PETN,  A  second  site  (labelled  2)  has  developed  in 
frame  8c.  By  the  final  frame,  fast  reaction  has  spread  throughout  most 
of  the  PETN . 

An  example  from  the  second  type  of  experiment  is  shown  in  Figure  9. 

In  this  case  a  single  polymer  disc  of  PS  was  pressed  with  the  PETN 
powder  to  form  a  composite  disc.  The  selected  frame  shows  ignition  at 
three  points  (I  to  3)  on  the  polymer/explosive  interface.  The  polymer 
is  the  lighter  central  area,  and  the  darker  annulus  is  the  PETN.  In 
the  lower  part  the  PETN  is  becoming  lighter  and  this  is  due  to  the  onset 
of  fusion.  At  its  outer  edge  the  PETN  is  beginning  to  produce  fast  jets. 
The  mottled  appearance  of  the  polymer  is  partly  due  to  catastrophic  failure 
and  partly  due  to  a  very  thin  layer  of  trapped  PETN.  This  is  why  reaction 
also  spreads  inwards  as  well  as  outwards.  Frame  9b  shows  a  later  state 
where  fast  reaction  is  more  developed.  The ^ expansion  velocity  of  the  PS 
disc  just  prior  to  ignition  wbb  ca.  35  m  s  .  The  catastrophic  failure 
occurred  between  frames.  Experiments  with  PC  showed  that  it  behaved 
in  a  very  similar  manner,  with  cracking  and  rapid  radial  expansion  followed 
by  ignition  at  the  interface. 

The  situation  with  PP  is  quite  different  and,  although  much  radial 
expansion  at  quite  high  velocities  is  observed,  the  polymer  does  not  crack 
or  shear  locally  and  ignition  does  not  occur. 

The  results  obtained  with  HMX/polymer  systems  are  similar  to  those 
obtained  with  PETN.  Figure  10  illustrates  an  impact  on  a  20  mg  sample 
containing  a  PP  disc  and  Figure  11  a  20  mg  sample  containing  a  PC  disc. 

In  both  cases  the  polymer  is  the  light  central  area  and  the  explosive  the 
darker  area.  Figure  10  shows  that  yield  and  flow  of  both  the  PP  and  HMX 
occur  but  no  ignition  results.  In  Figure  11  ignition  occurs  at  the  polymer/ 
HMX  interface  (frames  d  onwards)  after  catastrophic  failure  of  the  PC. 

The  findings  of  the  photographic  work  described  in  this  section  may 
be  summarised  as:  (i)  ignition  occurs  after  or  during  the  catastrophic 
failure,  cracking  and  fast  flow  of  the  polymer;  (ii)  the  initiation  site 
is  very  closely  associated  with  the  position  of  the  polymer  and  (iii)  of 
the  polymers  tested,  only  those  which  undergo  catastrophic  failure  are 
effective  in  promoting  ignition. 

2 . i  Impact  Behaviour  in  the  Drop-Weight  Test 

(a)  Impacts  of  polymer  samples 

Samples  of  the  polymeric  materials  were  impacted  in  a  drop-weight 
machine  at  strain  rates  similar  to  those  observed  in  the  photographic 
work  (t  ca.  I0J  s  *)so  that  differences  in  their  mechanical  behaviour 
could  be  observed.  The  drop-weight  machine  used  was  similar  to  that  des¬ 
cribed  by  Heavens  and  Field  (1974)  and  is  illustrated  schematically  in 
Figure  12.  The  sample  was  placed  between  rollers  Rl  and  R2  which  had  their 
edges  ground  down  to  give  a  sharply  defined  contact  area  of  diameter  5  mm. 
The  samples  used  were  discs  of  polymer  typically  1.5  mm  high  and  7  mm  dia¬ 
meter.  The  contact  area  was  lubricated  with  colloidal  graphite.  Two 
semiconductor  strain  gauges  were  mounted  in  series  diametrically  opposite 
to  each  other  on  R3  in  order  that  the  resistance  change  induced  in  the 
pair  of  gauges  during  impact  would  be  due  to  compressive  stresses  alone. 


The  system  was  calibrated  statically  using  an  Instron  machine. 


As  the  sample  diameter  was  greater  than  the  contact  area,  the  area  of 
contact  during  impact  remained  constant  throughout  the  whole  process.  The 
pressure-time  traces  were  obtained  on  paper-tape  utilising  a  Data  Lab  DL922 
transient  recorder  and  analysed  on  a  HP9825  calculator  using  a  programme 
based  on  the  method  of  analysis  of  Afanas'ev  and  Bobolev  (1971).  The 
programme  numerically  integrates  the  force-time  trace  produced  by  the 
.equipment  to  yield  a  velotiy-time  trace  for  the  weight 

V(t)  -  VQ-  M~ 1  C/0  F(t)  dt  (3) 

with  V0  the  velocity  of  the  weight  at  impact  and  M  its  mass.  A  second 
integration  gives  the  displacement  of  the  weight  as  a  function  of  time  z(t). 

z(t)  -  t/0  V(t)  dt  (4) 

z(t)  is  the  compression  of  the  whole  system  (sample  and  anvils)  and  a  small 
amount  i  (t)  due  to  the  compression  of  the  system,  which  can  be  determined 
by  performing  impact  with  no  sample  present,  is  subtracted  from  z(t)  to 
yield  the  compression  of  the  sample. 

z#(t)  ■  z(t)  -  zA(t)  (5) 

The  sample  strain  is  then  given  by 


e(t)  ■  In  -  (6) 

t  H  -  zfi(t)  i 

with  H  the  initial  height  of  the  sample  and  e(t)  the  strain  at'  time  t. 

The  stress  at  time  t  is  given  by  division  of  the  force-time  curve  by  the 
area  of  contact  A, 

o(t)  ■  F(t)/A 

Elimination  of  t  between  equations  5  and  6  yields  the  stress-strain  curve. 

Results  obtained  with  PC  and  PP  are  illustrated  in  Figure  13  (a)  and 
(b).  The  behaviour  of  PS  was  similar  to  that  of  PC.  The  main  difference 
between  the  two  materials  is  the  evidence  of  catastrophic  failure  in  PC 
(as  shown  by  the  sharp  drop  at  e  ■  1.8)  and  its  absence  in  PP.  Both  mater¬ 
ials  show  evidence  of  strain  softening  followed  by  orientation  hardening* 
but  the  strains  which  PP  is  capable  of  are  far  greater  than  those  achieved 
by  PC.  The  onset  of  orientation  hardening  also  takes  place  much  earlier  in 
PC  at  a  strain  of  ca.  0.5  as  opposed  to  the  value  of  2.0  observed  with  PP. 

At  strains  of  >2  the  thickness  of  the  PP  layer  becomes  comparable  to  the 
correction  factor  for  anvil  deflection  and  thus  the  calculated  strain 
become  increasingly  unreliable. 

The  shape  of  the  stress-strain  curve  may  therefore  be  inaccurate  at  high  strains 
These  results  confirm  the  findings  of  the  photographic  work  in  showing 
a  marked  difference  in  high  strain  rate  impact  behaviour  between  sensitising 
and  non-sensitising  materials  and  again  point  to  catastrophic  failure  as 
being  the  process  responsible  for  a  polymer's  sensitising  action. 
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(b)  Impact  behaviour  of  explosive /polymer  samples 

The  sensitivity  of  pure  explosives  PETN  and  HMX  and  explosive /polymer 
samples  were  compared  by  determining  their  50%  height  on  the  drop-weight 
machine  outlined  above.  The  30%  height  is  the  height  from  which  the  hammer , 
if  allowed  to  fall  freely,  would  have  a  30%  probability  of  causing  initiation. 
The  method  of  determination  was  based  on  the  "up  and  down"  method  which  is 
fully  described  in  Dixon  and  Massey  (1957). 

A  large  number  of  drop-weight  tests  were  carried  out  on  explosive/grit 
samples  to  determine  (i)  which  polymers  gave  rise  to  increased  sensitiveness , 
(ii)if  the  size  of  the  grit  particles  was  important ,(iii)  the  effect  of 
variation  of  the  percentage  weight  of  grit,  and  (iv)  if  a  polymer  which  sensi¬ 
tised  one  secondary  was  effective+in  sensitising  another.  All  the  experiments 
were  carried  out  using  twenty  25  -  1  mg  samples  of  powder  which  had  bean 
pressed  statically  to  10  kbar  to  give  pellets  of  5  mm  diameter  and  approxi¬ 
mately  0.8  mm  height.  In  these  experiments,  three  grades  of  grit  were  used 
which  have  been  called  coarse,  fine  and  very  fine.  In  the  'very  fins'  grade, 
all  the  particles  had  been  through  a  76  pm  mesh,  the  'fine'  grade  indicates 
that  the  grit  was  used  as  obtained  from  the  freezer-mill  with  a  maximum 
particle  mass  of  ca.  20  Ug.  The  coarse  material  was  as-received  grit  from 
AWRE  with  particles  of  mass  up  to  1  mg  and  up  to  I  mm  long.  These  reuslts 
are  set  out  in  Table  5.  Unless  otherwise  stated,  all  the  samples  contained 
20%  by  weight  of  grit 

Table  5:  30%  Heights  for  Explosives  and  Explosive/Binder  Samples 


Sample  50%  ht  (cm) 


PETN  (pure) 

21.8  ~  1.0 

PETN  +  PC  (very  fine) 

PETN  +  PC  (very  fine)  repeat 

10.5  ;  2.0  - 

10.7  -  0.6 

PETN  +  PS  (very  fine) 

PETN  +  PS  (fine) 

PETN  +  PS  (coarse) 

PETN  +  PS  (1  mm  high,  1mm  diam.) 

11.1  -  0.3 

13.0  i  1 .5 

10.4  *  0.6 

10.5  ±  3.0 

PETN  PMMA  (coarse) 

PETN  +  PMMA  (fine) 

21.0  t  1 .6 

19.2  t  0.8 

PETN  +  PP  (very  fine) 

PETN  +  PP  (very  fine)  10%  grit 

18.8  1  0.8 

19.2  t  0.7 

HMX  (pure) 

29.0  t  2.0 

HMX  +  PC  (fine) 

16.0  -  1.5 

HMX  +  PS  (very  fine) 

15.1  i  1.0 

HMX  +  PMMA  (coarse) 

22.6  -  1.0 

HMX  +  PP  (fine) 

24.4  -  0.8 

AgN^  (pure) 

19.3  -  0.8 

AgNs  +  PC  (fine) 

16.8  ±  1.3 
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From  the  results  set  out,  it  can  be  seen  that  the  size  of  the  grits 
used  does  not  have  any  great  effect  on  the  results  produced.  It  might  be 
expected  that,  if  the  effectiveness  of  the  grit  demanded  on  its  ability  to 
disturb  the  flow  uf  the  explosive  as  it  was  extruded  from  between  the 
anvils,  the  lnrgi  particles  by  virtue  of  their  size  would  prove  more  effec¬ 
tive,  the  very  fine  particles  being  so  small  as  to  be  carried  along  in  the 
flow  of  the  explosive.  The  results  obtained  show  that  this  is  not  a  major 
factor. 

The  results  also  show  that  the  grits  are  not  specific  in  their  action: 

1'C  and  PS  being  equally  effective  in  sensitising  HMX  and  '’ETN,  while  PP 
and  PMMA  only  slightly  sensitise  both  of  these  secondary  explosives.  As 
noted  before  (see  for  example,  Heavens,  1973,  and  Heavens  and  Field,  1974) 
primary  explosives  such  as  lead  azide  PbCN*)*  and  silver  azide  AgN^  are  not 
particularly  sensitive  in  the  drop-weight  test  provided  the  anvils  are  well- 
aligned.  However,  as  the  table  shows,  the  polymers  do  have  a  significant 
effect  on  the  sensitivt ness  of  AgN.  and  in  a  consistent  way.  Finally,  in 
one  set  of  experiments  the  percentage  of  grit  was  reduced  to  10Z  but  this 
caused  no  significant  change. 

2 • "  Thermal  Properties 

(a)  Mechanical  effects  and  thermal  properties 

The  mechanical  failure  properties  of  polymer  additives  appear  to  be 
of  major  importance  in  determining  their  sensitising  effect.  Those  mater¬ 
ials  which  undergo  bull  plastic  deformation  proving  ineffective  as  sensitises 
and  those  which  fail  ct  tastrophically  and  with  intense  local  deformation 
effective.  Thermal  pre ierties  areof  major  significance  in  determining 
the  failure  modes  of  pc Lymers  particularly  when,  high  strain  rates  are  in¬ 
volved.  This  is  because  they  soften  at  relatively  low  temperatures  and  have 
low  densities,  specific  heats  and  thermal  conductivities.  Additionally, 
materials  can  dcsensiti?e  explosives  by  effectively  acting  as  heat  sinks 
and  quenching  hot  spots  before  they  can  fully  develop,  The  work  of  Bowers, 
Romans  and  Zisman  (1973)  for  example,  indicates  that  sensitivity  decreases 
with  increasing  specific  heat  of  additive,  those  additives  which  are  most 
effective  in  this  respect  being  materials  such  as'Superla  wax"or "Carbowax" 
which  have  low  melting  points  and  thus  absorb  their  heat  of  fusion. 

In  cases  where  thermal  properties  were  not  known  they  were  measured 
using  a  Perkin-Elmer  DSC-2  calorimeter. 

The  results  obtained  show  that  both  PC  and  PS  have  specific  heats 
which  range  from  1.25  kJ  kg-1  K-'  at  room  temperature  to  2.35  kJ  kg-1  K“’ 
at  600  K,  the  increase  jeing  reasonably  linear.  PC  has  a  melting  point 

of  420  to  430  K  which  oily  introduces  a  very  slight  increase  in^he  effective 

specific  heat,  giving  a  latent  heat  of  fusion  of  ca.  7.2  k.I  kg  .  PS  shows 
melting  over  the  range  4 70—500  K  with  a  heat  of  fusion  of  ca.  6,6  kJ  kg  '. 

PP  has  a  greater  specific  heat,  varying  from  1.92  kJ  kg”'  K”'  at  room  temp¬ 

erature  to  2.09  kJ  kg-'  K-'  at  450  K.  It:  melted  over  the  range  390  - 
460  K  with  a  latent  hea:  ot  fusion  of  ca.  64.0  kJ  ka“'  K“' .  PMMA  had  a 
relatively  constant  spe:ific  heat  of  2.9  kJ  kg”'  K“*  over  the  range  470  - 
630  K  where  a  large  dec  imposition  endotherm  was  produced.  From  the  results 
set  out  above,  it  could  be  predicted  that,  on  the  basis  of  thermal  properties 
alone,  the  order  for  ho:-spot  quenching  for  the  polymers  investigated  would 
be  PP,  PMMA,  PC,  PS. 
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Measurements  were  also  made  on  the  specific  heat  and  latent  heat  of 
fusion  of  PETN.  It  was  found  that  specific  heat  remained  almost  constant 
at  1.09  kJ  kg  up  to  the  melting  temperature  of  414  K.  The  latent  heat 
was  determined  to  be  144.7  kJ  kg-*  with  decomposition  occurring  rapidly 
after  melting.  A  summary  of  the  thermal  properties  of  the  polymers  and  ex¬ 
plosives  used  in  this  work  is  given  in  Table  6  . 

2.b  Polymer  Friction  Experiments 

The  purpose  of  the  experiments  described  in  this  section  was  to  measure 
the  temperatures  that  a  polymer  reaches  during  deformation.  In  order  to  do 
this,  a  friction  machine  was  constructed  which  allows  the  rubbing  polymer 
interface  to  be  viewed  by  an  infra-red  detector.  Friction  cannot  be  regarded 
as  just  a  surface  shearing  effect  since,  if  one  of  the  materials  is  harder 
than  the  other,  the  asperities  on  the  harder  surface  can  plough  out  the  softer 
material  to  an  appreciable  depth  below  the  surface.  It  is  thus  possible  for 
considerably  and  varied  deformations  to  take  place  in  a  frictional  situation, 
and  estimates  of  the  surface  temperatures  achieved  are  of  interest  both  to 
the  present  investigation  and  the  subject  of  polymer  friction  and  wear  in 
general . 

(a)  Apparatus  and  technique 

The  apparatus  constructed  to  perform  these  experiments  was  basically 
a  modification  of  that  described  by  Bowden  and  Thomas  (1954)  and  is  shown 
schematically  in  Figure  14.  It  consisted  of  a  50  mm  diameter  rotating  sapphire 
disc  (infra-red  transmitting  .to  ca.  5.3  pm)  against  which  a  polymer  pin  was 
loaded  by  means  of  a  balance  arm.  The  area  of  the  interface  under  observation 
was  defined  by  a  0.8  mm  slit  mounted  directly  above  the  rotating  disc.  The 
disc  rotated  at  a  frequency  of  ca.  56  Hz  giving  an  interfacial  velocity  of 
ca.  7ms*"1.  Above  this  disc  was  a  PMMA  chopper  with  radial  regions  in  the 
sequence} blacked-out,  slit,  PMMA,  slit,  blacked  out  etc.  The  disc  was  ro¬ 
tated  so  that  the  radiation  from  the  interface  was  chopped  at  a  frequency 
of  10  kHz.  This  allowed  temperature  variations  of  duration  1  ms  or  greater 
to  be  followed  (since,  for  any  one  temperature  determination,  a  complete 
blacked-out  to  blacked-out  sequence  must  be  obtained,  i.e.  5  ’chops').  After 
chopping,  the  radiation  was  detected  by  a  Mullard  RPY  36  infra-red  detector 
which  is  sensitive  over  the  range  1  -  5.5  pm  with  a  sharp  cut-off  at  5.5  pm. 
The  detector  was  biased  to  give  the  best  signal-to-noise  ratio.  The  detector 
output  was  taken  to  a  Tektronix  7A22  amplifier  and  the  amplifier  output  was 
fed  to  a  Datalab  DL922  transient  recorder.  Traces  were  then  punched  onto 
paper-tape  for  analysis  on  the  Hewlett-Packard  HP  9825  desk-top  calculator 
and  plotter. 

The  purpose  of  the  PMMA  filters  was  to  provide  an  absorption  band 
(ca.  3.5  pm)  which  would  absorb  a  differing  proportion  of  the  signal  depend¬ 
ing  on  the  temperature  of  the  source.  Thus,  by  calculating  the  ratios  I/Ip 
(with  I  the  intensity  through  the  slit  and  Ip  through  the  filter)  for  sources 
of  various  known  temperatures,  a  calibration  curve  of  I/Ip  against  temperature 
could  be  drawn  and  estimates  of  temperatures  of  other  sources  made  using  the 
curve.  In  order  for  the  calibration  to  be  valid,  it  was  assumed  that  the 
sources  (calibration  and  specimen)  were  behaving- as  black-bodies.  Problems 
associated  with  different  substances  having  differing  black-body  emissivities 
arc  eliminated  by  the  use  of  the  ratio  I/Ip. 
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(b)  Calibration 

The  theory  of  calibration  follows  the  treatment  of  Land  ( 1944)  and 
Parker  and  Marshall  (1948).  From  Planck's  radiation  law  we  have 

B,  _  -  C  e,  X"5  /  (eB/XT-  1)  (7) 

a,  r  a 

with  EXjt  ^e  energy  emitted  by  the  black-body  at  wavelength  \  and  temperature 
T,  C  and  B  constants  and  ex  the  eraissivity  of  the  surface  at  wavelength  A. 
This  can  be  approximated  to  Wien's  formula. 

E  •=  A  G  ex  t  A-5  exp(-B/A  T)  (8) 

P  P  P 

with  E  the  radiant  energy  received  by  the  detector , A  a  constant  dependent 
on  the  area  of  the  source,  G  a  constant  dependent  on  the  geometry  of  the 
system,  t  the  transmission  coefficient  of  the  sapphire  and  X  the  effective 
wavelength  of  response  of  the  detector.  The  relation  between  the  response 
of  the  detector  H  and  the  energy  falling  on  it  is  of  the  form 


with  a  a  constant. 


Thus  from  equations  8  and  9: 

-5 

In  H  *  m(ln  (A  G  t)  +  in(e.  A  )  -  B/X  I)  (10) 

p  p  p 

When  the  radiation  passes  through  thecPMMA  filter  the  optical  properties 
of  the  system  are  changed  (X,F.,A  and  Xpchange)  but  the  factors  A,  G  and 
t  remain  the  same.  Thus  if  H'  p  is  the  reduced  output  on  passing  through 
•the  filter  we  may  write 

In  H*  ■  m(ln(AGt)  +  ln(e’  A’"5)  -  B/X'T)  (11) 

P  P  P 

From  equations  10  and  11 


02) 


This  relation  is  independent  of  the  source  area  and  geometrical  fractors  and  is 
independent  of  the  absolute  values  of  the  emissivities.  It  depends  only  on 
the  differences  in  optical  properties  of  A  and  -A'  (and  this  is  dependent  only 
on  the  properties  of  the  system  with  and  w?thout  the  filter).  Thus  a  plot 
of  In-,  against.  J.  should  give  a  straight  line  and,  in  fact,  this  is  what 
is  fotind .  T 


Initially  calibration  was  carried  out  using  a  2.5  mm  thick  PMMA  chopper 
disc  and  a  0.27  mm  diameter  platinum  wire  as  an  artificial  hot-spot.  The 
temperature  of  the  wire  was  varied  by  changing  the  DC  current  flowing  through 
it.  The  voltage  across  a  central  portion  of  the  wire  was  monitored  via  two 
25.4  pm  diameter  platinum  wires  which  were  spot-welded  to  the  larger  wire  48.0  mm 
apart.  These  wires  were  thin  enough  not  to  affect  the  temperature  of  the  larger 
wire  which  was  considered  to  be  uniform  over  the  48.0  mra  long  middle  portion. 

From  the.  knowledge  of  the  resistance  of  this  central  portion,  its  temperature 
was  determined  using  standard  tables.  The  central  portion  was 
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(c)  Results 

In  the  experimental  investigation ,  polymer  pins  of  two  different  dia¬ 
meters  were  used:  small  pins  of  diameter  ea.  1.5  nun  and  larger  pins  °f  dia¬ 
meter  ca.  6  mm.  The  smaller  diameter  pins  were  used  in  the.  initial  experiments 
but  for  the  later  measurements  large  pins  were  used  since  the  smaller  pins 
gave  a  poor  signal-to-signal  noise  ratio  and  were  very  easily  deformed  due 
to  their  large  length/diameter  ratio.  The  main  problem  with  the  use  of  the 
larger  pins  was  that  considerable  heating  of  the  sapphire  disc  took  place 
and  surface  temperatures  of  up  to  360  K.  were  recorded  oh  placing  a  thermo¬ 
couple  in  contact  with  the  disc  immediately  (1  -  2  s)  after  doing  a  run. 

Due  to  the  considerable  thermal  mass  of  the  disc.,  its  rate  of  cooling  was 
quite  slow  (10  -  20  degrees  per  minute)  mid  so  the  figure  quoted  above  gives 
a  good  indication  of  the  bulk  temperature  reached  by  the  sapphire,  in  the 
course  of  an  experiment. 

An  experiment  in  which  the  signal  level  was  recorded  with  the  large  pin 
in  contact  with  the  disc  and  again  immediately  after  its  removal  indicated  that 
the  energy  due  to  the  bulk  heating  of  the  disc  accounted  for  approximately 
one  third  of  the  total  signal.  Such  a  large  contribution  from  a  lower 
temperature  source  will  give  a  mixed  signal  which  will  tend  to  reduce  the 
measured  black-body  temperature  to  a  value  near  the  mean.  Thus  the  readings 
obtained  with  the  larger  pins  will  probably  be  underestimates  of  the  interface 
temperature.  In  the.  table  of  results  given  below  (Table  V)  the  maximum  black- 
body  temperature  recorded  by  the  system  for  any  given  polymer,  irrespective  of 
the  type  of  pin  used  is  quoted. 

The  results  will ,  of  course  »  be  averages  over  the  total  portion  of  the 
interface  area  within  the  detector's  field  of  view  and  it  is  quite  likely  that 
some  regions  (those  being  highly  deformed)  will  have  temperatures  above  this 
average  value  and  others  will  have  lower  values.  By  far  the  largest  amount 
of  experimental  time  was  spent  on  the  polymers  ?P  and  PC  so  that  with  an 
increased  amount  of  observation  time  it  may  be  possible  to  record  higher 
maximum  values  of  temperature  for  the  three  other  materials.  Ic'  is,  however,  in¬ 
teresting  that  the  temperatures  obtained  follow  the  order  obtained  for  the 
sensitiveness  of  the  polymer/explosive  mixtures.  It  is  also  clear,  particularly 
in  the  case  of  PC,  that  temperatures  well  above  the  commercially  quoted 
softening  temperature  of  4 10  K  and  the  observed  'melting'  temperature  (in 
previous  work)  of  430  K  are  obtainable.  In  the  case  of  PF  however  this  was  not 
found  to  be  the  case  and  no  temperatures  above  the  commencement  of  melting 
at  390  K  were  obtained. 

Table  7:  Temperatures  Achieved  During  Polymer  Friction 


Polymer _ _ Maximum  Interfacial  Tcmperature/K 


pp 

360  -  30 

PMMA 

420  *  30 

PS 

450  -  20 

PC 

490  -  10 

(c)  Results 


i‘>. 


In  lilt’  experimental  inves  li  g.iL  ton ,  polymer  pins  of  !  wo  different  dia¬ 
meters  were  used:  small  pins  of  d  iuinet.-r  ra.  I  .  mm  and  larger  pins  of  dia- 
motor  <  a .  6  mm.  Tin  smaller  diameter  pint;  worn  used  in  the  ini  tial  oxpor  i merit 
but  for  the  later  moa  summon  i  s  large  pins  woe  used  since  the*  smaller  pins 
gave  a  poor  signal -t ('-signal  noise  ratio  and  were  very  easily  deformed  due 
to  their  large  lrngth/d Lame ter  ratio.  The  main  problem  with  the  use  of  the 
larger  pins  was  that  considerable  boat  ing  of  the  sapphire  disc  took  place, 
and  surface  temperatures  of  up  to  360  K  were  recorded  on  placing  a  thermo¬ 
couple  in  contact  with  the  disc  immediately  0  -  ?.  s)  after  doing  ft  run. 

Due  to  the  considerable  thermal  mass  of  the  dine,  its  rate  of  cooling  was 
quite  slow  (10  -  20  degrees  per  minute)  and  so  the  figure  quoted  above  gives 
a  good  indication  of  the  bulk  temperature  reached  hv  the  sapphire  in  the 
course  of  an  experiment  . 

An  experiment'  in  which  the  signal  level  was  recorded  with  the  large  pin 
in  contact,  with  the  disc  and  again  immediately  after  itfl  removal  indicated  that 
the  energy  due  to  the  bulk  heating  of  the  disc  accounted  for  approximately 
one  third  of  the  total  signal.  Such  a  large  contribution  from  a  lower 
temperature  source  will  give  a  mixed  signal  which  will  tend  to  reduce  the 
measured  black- body  temperature  to  a  value  near  the  mean.  Thus  the  readings 
obtained  with  the  larger  pins  will  probably  he  underestimate*  of  the  Interface 
temperature.  In  the  table  of  results  given  be*. 'Jow  (Tabic:  7)  the  maximum  hlack- 
body  temperature  recorded  by  the  system  for  any  given  polymer,  irrespective  of 
t.'he  type  of  pin  used  is  quoted. 

The  results  will,  nf  course  >  be.*  averages  over  the  total  portion  of  the 
interface  area  within  the  detector's  field  oi.  view  and  it  is  quite  likely  that 
Koine  regions  (those  being  highly  deformed)  will  have  tempo  rot  urea  above  this 
average  value  mid  others  will  have  lower  values,  by  far  the  largest,  amount 
of  experimental  time  was  spent  on  the  polymers  FT  and  PC  ro  that  with  an 
increased  amount  of  observation  time  it  may  be  possible'  to  record  higher 
maximum  values  of  temperature  for  the  three  other  material s ,  IV  is,  however,  in¬ 
teresting  that  the  temperatures  obtained  follow  the  order  obtained  for  the 
seriNitlvoncsJfl  of  the  polymer /cxplufii vo  mixtures.  Tt  is  also  clear,  particularly 
in  the  case  of  PC,  that  tempera t .urea  well  above  the  commercially  quoted 
softening  temperature  of  610  K  and  the  observed  'melting'  temperature  (in 
previous  work)  of  430  K  are  obtainable.  In  the  case  of  PI’  however  this  war.  not 
found  to  bo  the  case  and  no  temperatures  above  the  commencement  of  melting 
at  IDO  K  were  obtained. 


Table  /:  Temperatures  Achieved  During  P£b}JJ!JJ?r  Pri  ct^ion 

Maximum  Inter facial  Temper at ure/K 


pp 

360  *  JO 

PMMA 

420  *  30 

PS 

4  r.p  -  :>o 

PC 

4 SO  -  10 

Examination  of  the  polymer  pins  after  the  frictional  experiments  reveal 
that  in  the  case  of  PC  and  PS  severe  deformation  of  the  polymer  had  been 
produced  with  evidence  of  much  stretching  and  shearing  of  molecular  layers. 

The  PP  pin,  however,  present  a  much  more  rounded  appearance  which  indicated 
that  melting  of  the  surface  layers  hod  taken  place. 

Figure  17  shows  the  overall  appearance  of  the  pins  after  the  experiments 
(a)  is  a  PC  pin  and  the  many  layers  which  have  sheared  off  and  built  up  on 
the  trailing  edge  are  clearly  visible,  (b)  is  a  PP  pin  and  shows  the  rounded, 
meltcd-looking  appearance  of  the  surface,  (c)  is  a  PMMA  pin  and  shows  the 
roughened  surface  obtained  with  no  build  up  on  the  trailing  edge.  All  the 
pins  in  Figure  17  were  of  diameter  6  mm. 

Microscopic  examination  confirmed  that  very  severe  deformation  of  PC  had 
taken  place  with  clear  evidence  of  considerable  strain  in  the  frictional 
direction  (Figure  18(a)).  PMMA  (Figure  18(b))  presented  a  rather  'powdery' 
appearance  with  many  particles  of  material  on  the  surface.  The  wearing  away 
of  the  PMMA  pins  proceeded  in  a  more  brittle  manner  with  small  particles  chipped 
out  of  the  surface  and  with  no .accumulation  of  sheared  or  melted  layers  on  the  . 
trailing  edge  of  the  pin.  The  surface  of  PP  on  the  other  hand  was  relatively 
featureless;  as  might  be  expected  from  a  melted  and  refrozen  layer. 

Although  there  is  clear  evidence  that  PC  had  deformed  by  the  shearing  off 
of  layers,  there  was  no  evidence  of  catastrophic  failure  having  taken  pla  i: 
the  shearing  could  be  followed  by  the  naked  eye  and  took  place  over  a  time- 
scale  of  ca.  1  s  for  a  6  mm  pin.  Nonetheless,  temperatures  well  above  the 
commerical  softening  temperature  were  recorded  with  PC,  and,  when  allowance 
is  made  for  the  fact  that  these  were  averages  over  the  whole  surface,  it  does 
seem  possible  that,  under  conditions  of  extreme  (and  probably  locally  rapid) 
deformation,  very  high  temperatures  could  be  achieved  by  polymers  such  as  PC. 
These  temperatures  would  be  well  above  the  limits  set  by  the  old  ideas  of  the. 
temperature  rise  achieved  in  friction  being  limited  by  the  melting  point  of 
the  lower  melting  point  material  which  was  based  on  the  frictional  properties 
of  metals.  Any  increase  in  melting  point  of  the  polymers  due  to  the  force 
applied  by. the  counterweight  would  be  negligible  since  the  maximum  weight 
used  with  the  6  mm  pins  was  2  kg  and  with  the  1.5  mm  pins  was  500  g  giving  a 
maximum  stress  of  3  MPa,  which  corresponds  to  ca.  30  bar.  This  would  have  an 
insignificant  effect  on  the  melting  point,  variations  of  ca.  200  bar  being 
needed  to  cause  an  appreciable  (ca.  10°  K)  change  (Zoller  1978).  It  does 
therefore  appear  that  high  temperatures  can  be  reached  during  the  deformation 
of  polymers,  and,  in  the  design  of  any  system  containing  a  polymer,  consideration 
must  be  given  to  the  possibility  of  high  temperature  regions  during  deformation. 


2.7  Conclusions  and  Discussion 


A  combinati' 
thermal  analysis 
action  of  cxplos 
tested,  were  foui 
emiicmn  in  ini  pa 
after  rapid  and  : 
presence  of  sens 
necessary  fur  ip. 
ndtli.  t  ivo  are  imp 
having  a  comb  in* 
of  fusion,  wine1 


in  of  high-speed  photographic,  pressure  measuring  and 
techniques  have  been  used  to  study  the  behaviour  and  inter¬ 
nes  and  polymers.  Chemical  interactions,  for  the  materials 
'd  to  be  of  minor  importance  compared  with  mechanical  consid- 
t  sennitisation.  Pure  explosives  were  found  to  ignite  only 
evere  deformation  of  the  sample  had  occurred.  In  the 
tising  grits  however,  this  severe  deformation  war,  not 
it  ion.  Both  f.lie  .fracture  and  thermal  properties  of  an 
t  tan  L  in  del  ermi  ni.ng  its  sensitivity  effect.  A  polymer 
inn  of  high  si renglh  with  low  specific  heat  and  low  heat 
lend  to  promote  catastrophic  failure, will  be  likely  to 


cause  sensitisation.  In  order  to  teat  this  hypothesis,  two  polymers  were 
chosen  on  the  basis  of  their  strength. and  thermal  capacity;  polyester  which 
has  high  strength  and  low, capacity  and  high  density  polythene (HOPE)  which 
has  a  low  strength  and  high  capacity.  They  were  mixed  with  PETN  and  tested 
in  the  drop-weight  machine  in  , the  same  manner  as  the  samples  described  earl  ier. 
The  results  obtained,  together  with  a  few  comparison  results  from  Table  5  are 
set  out  in  Table  8.  Thermal  properties  can  be  found  in  Table  6. 


Table  8:,  50%  Heights  for 


nor  Samples 


Sample  ___ _ . . _  •;  50%  height  (cm) 


PETN  (pure 

21.8  - 

1.0 

PETN  +  HDPE 

24.8  - 

1  .5 

'  i  '  ■ 

PETN  ♦  PP 

18.8  - 

0.8 

PETN  +  Polyester 

15.5  * 

1.0 

• 

PETN  +  PS  :  !IJ  -  0.3 


These  results  confirm  the  predictions  that  polyester  should  sensitise 

and  HOPE  should  not.  In  fact  it  appears  that  HOPE  is  acting-  as'  a  desehsitiscr. 

I  .  ,  I  ...  • ;  • 

The  conclusion  reached,  therefore,  is  that  polymers  which  are  most 
likely  to  sensitise  are,  those  such  as  PC  which  are  tough  below  their  glass 
transition  temperature  (as  opposed  to  PMMA  which  behaves  in  a  much  more 
brittle  manner),  and  are  likely  to  fail  catastrophically.  Those -least 
likely  to  sensitise  are  those  which  have  , a  low  glass  transition  temperature 
(PP  or  HDPE)  and  which  deform  in  bulk  with  little  local  concentration  of 
energy.  Polyester  is  somewhat  intermediate  in  behaviour  (and  sensitising' 
effect)  between  PC  and  PMMA  as  it  it  tougher  than  PMMA  but  much  more  brittle 
than  PC,  and  so  fits  quite  well  into  the  expected  pattern. 


In  the  case  of  the  former  the  melting  point  of  the  particle  is  the 
important  property  (Bowden  and  Gur ton  1949).  That  this  is  so  is  essen¬ 
tially  because  they  are  materials  with  well-defined  melting  points  and 
relatively  large  latent  heats;  the  hot  spot  is  generated  to  the  melting  . 
point.  In  contrast  the  polymeric  materials  which  sensitise,  soften  over 
a  range,  of  temperatures  anc!  have  low  latent  heats  and  thermal  conductivities. 
If  heat  is  produced  by  localised  catastrophic1  failure  hot  spots  will  form. 
However,  when  softening  of  the  polymer  starts  the  material  will  still 
maintain  a  high  viscosity.  Continuing  deformation  will  be  able  t:o  con¬ 
tribute  further  temperature  rises  by  viscous  heating.  Thor-  thermal 
losses  in  the  case  of  a  polymer  arc  very  low  compared  with  most  other 
materials  because  of  the  low  latent,  heat,  thermal  capacity  and  conductivity, 
and  the  hot  spot  is  not  quenched  so  readily. 

The  friction  experiment  gave  evidence  of  temperatures  well  above  the 
softening  point.  In  this  experiment,  an  average  temperature  for  a  contact 
area  of  ca.  25  mm2  was  maintained.  Local  hot  spot  temperature  would  be 
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expected  to  exceed  the  valued  recorded.  In  other  experiments  in  the  laboratory, 
the  temperature  rise  at  the  tip  of  a  crack  propagating  in  a  polymer  was  re¬ 
corded,  again  using  I  R  techniques  (Fuller,  Fox  .and  Field  l97!i).  Two  impor¬ 
tant  results  emerged  from  this  work.  Tlu  first  was  that  temperatures  wall 
above  the  softening  points  were  observed.  Secondly, that  the  fracture  surface 
energy  increased  with  crack  velocity.  This  latter  re.au]  1.  was 

explained  by  the  high  temperatures  causing  softening  of  the  material  at  the 
crack  tip  thus  requiring  more  energy  to  bo  expended  in  crack  growth. 

It  is  clear  from  the  earlier  work  and  the  present  work  with  explosive/ 
polymer  samples  that  some  polymerB  can  reach  very  high  local  temperatures 
during  their  deformation.  This  is  clearly  relevant  to  the  safe  handling 
of  explosives  since  polymeric  muterials  may  be  used  for  packaging  or  ns 
binding  materials  in  explosive  composites.  The  results  alHo  have  interest 
to  studies  of  the  fracture,  friction,  erosion  and  wear  of  polymers. 

ACKNOWLEDGEMENTS 


This  work  was  sponsored  in  part  by  AWRE,  Aldermaston  and  in  part  by 
the  US  Government  through  its  European  Research  Office.  Grants  from 
both  sponsors  are  gratefully  acknowledged. 

REFERENCES 


Afauas'ev  G.T.  and  Bobolav,V.K. ,  1961,  Dolk.  Akad.  Nauk.  S'SSR  138,886. 

Afanas'ev  G.T. ,  Bobolev,  V.K.  Kazassova,Yu  A  and  Knrabam>v,Yu  F.,  1972, 

Fiz,  Gorenia  i  Vzryva  8(2)  229. 

Afanas'ev  G.T.  and  Bobolev  V.K.,  1971,  Israel  Program  for  Scientific 
Translations  "Initiation  of  solid  explosives  by  impact" 

Andreev,  K.K.,  Maurina,  N.D.  and  Rusakovu,  Yu.  A.,  1955,  Dokl..  Akad.  Nauk, 

SSSR,  J05,  533. 

Bean,  C.M.,  1973,  Private  Communication,  AWRE. 

Bolkhovitinov,  L.G.,  Dokl.  Akad.  Nauk.  SSSR,  \25_,  570, 

Bolkhovitinov,  L.G.  and  Pokhil,  P.F.,  1958,  Aok)  .  Akad.  Nauk.  SSSR,  123,  637. 

Bowden,  F.P.,  Proc.  Roy.  Soc.  1950,  204 

Bowden,  F.P.,  1963,  9th  Int.  Symp.  Academic  Press,  N.  York. 

Bowden,  F.P.  and  Yoffe,  A.D.,  1949,  3rd  Symp.  on  Combustion  and  Flame  Phenomena, 
p.  551,  Baltimore!  Williams  and  Wilkins  Co. 

Bowers,  R.C.,  Romans,  J.B,  and  Eisman,  V.A.,  1973,  I.  and  F.C  Report,  12,2, 

Bowden,  F.P.  and  Gurton,  O.A.  ,  1949,  Proc.,  Roy.  Soc.  A  I9B,  337. 

Bowdui,  F.P,,  Mulcahy,  M.F.R.,  Vines,  R.G.  and  Yoffe,  A.D.,  1947,  Proc. 

Roy.  Soc.  A_188,  291. 

Bowden,  F.P.  and  Thomas,  P.lt.,  1954,  Proc.  Roy.  Soc.  A223,  29. 

Bowden,  F.P,  and  Yoffe,  A.D.,  1952,  CUP,  "Initiation  and  growth  of  explosion 
in  liquids  and  solids". 

Bowden,  F.P.  and  Yoffe,  A. I).,  1958,  But tnrworths /'Fast  reaction  in  solids". 
Chnudhri ,  M.M.  and  Field,  J.E.,  1974,  Proc.  Roy.  Soc.  A340,  113. 


23. 


Cundall ,  R.B.,  Palmer,  T.F.  and  Wood,  C.E.,  1978,  J.  Chem.  Soe.  Faraday 
Trana.,  I  74(6),  1339. 

Dixon,  W.J,  and  Massey  F.J.,  1957,  McGraw-Hill  "Introduction  to  statistical 
analysis". 

Eirich,  F.R.  and  Tabor,  D. ,  1947,  Proc.  Camb.  Phil.  Soc.  4|4,  566. 

Fox,  P.G.,  1970,  J.  Solid  Stats  Chen.  2,  491. 

Freeman,  E.S.  and  Carroll,  B.,  1958,  J.  PhyB.  Chem,  6_2 ,  394. 

Fuller,  K.N.C.,  Fox  P.G.  and  Field,  J.E.,  1975,  Proc..  Roy.  Soc.  A341 ,537. 
Gurrieri,  S.,  Siracusa,  G.  and  Cali,  R. ,  1974,  J.  Thermal  Analysis  6_,  293. 
Hauser,  H.M.  and  Field,  J.E.,  1978,  Thermochimica  Acta  2£,  I. 

Heavens,  S.N.,  1973,  University  of  Cambridge,  Ph.D.  Thesis. 

Heavens,  S.N.  and  Field,  J.E.,  1974,  Proc.  Roy.  Soc.  A3 38,  77. 

Kholevo,  N.A. ,  1946,  Tr.  Kazansk  Khim.  Tech.  Inst.  _[0,  91. 

Land,  T.,  1944,  J.  Iron  and  Steel  Inst.  |49,  481. 

Maycoek,  J.N.  and  Vernecker,  V.R,  Pai,  1970,  Thermochimica  Acta  J^,  191. 

Ng.  W.L.,  Field,  J.E.  and  Hauser,  H.M.,  1976,  J.  Cham.  Soc.  Perkin.  Trans. 
II,  637, 

Porker,  R.C.  and  Marshall,  P.R.,  1948,  Proc.  Inst.  Mech.  Engrs,  158,  209. 

Kanby,  B.G.  and  Uabek,  J.F.,  1975,  Wiley,  "Photodegradation,  Photo-oxidation 
and  Photodestabiliaation  of  Polymers" 

Reich,  L. ,  1973,  Thermochimica  Acta,  £,  433. 

Rideal,  E.K.  and  Robertson,  A.J.B.,  1948,  Proc.  Roy.  Soc.  AI95,  1 35 
Robertson,  A.J.B.,  1948,  .J.  Soc.  Chem.  Ind.  67,  221. 

Rogers,  R.N.  and  Morris,  E.D.,  1966,  Analyt.  Chem.  38,  412. 

Taylor,  W.  and  Weale,  A.,  1932,  Proc.  Roy,  Soc.  A.  1 38,  92. 

Ubbelohde,  A.R.,  1948,  Phil.  Trans.  A241,  199.' 

Walker,  F.E.  and  Croon  Le  Roy,  1976,  US  Government  Report,  AD-A022812 
Winter,  R.E.  and  Field,  J.E.,  1975,  Proc.  Roy.  Soc.  A343,  399. 

Zoller,  P.,  1979,  J.  Appl.  Polymer  Sci.  23,  1051. 


i  x  ■ Tx~.~- :.TFX -TrT.TrrTrT  :■ ■ 


24. 


FIGURE  CAPTIONS 

1.  Activation  energy  plot  for  the  thermal  decomposition  of  PETN  for 
various  values  of  reaction  order  n.  The  n  value  which  gives  the  best 
straight  line  segments  is  one  appropriate  to  the  reaction  (here  0.8). 

These  two  steps  (see  dashed  line)  corresponding  to  sublimation,  stage  T, 
and  decomposition  stage  11. 

2.  Experimental  arrangement  at  instant  of  impact,  W;  drop-weight;  G, 
glass  blocks;  M,  mirrors;  S,  sample.  The  upper  glass  block  is 
attached  to  the  weight. 

3.  Impact  on  a  disc  of  PP.  The  sample  deforms  by  bulk  plastic  flow 
throughout.  Diameter  of  fiel-d  of  view  20  mm,  Frame  time/ps  : 

a, 0;  b,70;c,266.  The  zero  times  given  in  this  and  subsequent  sequences 
do  not  correspond  to  time  after  initial  contact. 

4.  Impact  of  a  disc  of  PS.  The  sample  deforms  initially  by  bulk  plastic 
flow  (frames  a  and  b)  but  then  deforms  catastrophically  (frame  c). 

Diameter  of  field  of  view,  2,!i  run.  Frame  times/ps  :  a,0;b,2l  ;c,28. 

5.  Diameter  versus  time  traces  for  impacts  on  PP  and  PS.  Note  that 

the  diameter  of  the  l'S  has  approximately  doubled  before  the  entastropic 
failure,  but  that  the  PP  disc  has  a  final  diameter  several  times  its 
initial  value, 

6.  Micrographs  of  a  deformed  PS  disc  showing  that  the  catastrophic  failure 
is  a  result  of  fracture  and  localised  shear.  The  enlarged  view  in  (b) 
in  from  region  S  in  (a).  The  view  has  been  rotated  so  that  the  shear 
hands  are  horizontal. 

7.  Impact  on  a  pellet  of  PETN.  Frame  a  allows  the  5  mm  diameter  pellet 
before  Impact.  Rapid  flow  mid  jetting  occur  as  for  powdered  layers. 

Fusion  of  the.  layer  starts  in  frame  d.  The  whole  layer  is  transparent 
und  flowing  rapidly  hy  frame  g  when  the  first  to  two  ignition  sites 
forms,  The  circular  feature  in  frame  h  is  an  artifact.  Diameter  of 
field  of  view  20  mm.  Frame  times/ps  :  a ,0;b,67;c ,22 1 ;d,308;e, 389 ; 

f  ,469jg,47b;h,482;i489. 

8.  PETN  pellet  surrounded  by  5  PS  discs  Sec  text  for  details,  Two  ig¬ 
nition  sites  at  the  polymer/explosive  interface  are  labelled  in  frames 
d  and  c.  Diameter  of  field  of  view  20  nun.  Frame  times/us  :a,0; 

b ,  9 1 ; c 1 40 ; c , 1 6 1 ; c , I b8 ; f ,  I  75 . 

9.  25  mg  disc  of  PETN  containing  a  PS  disc.  Three  ignition  sites  are 
labelled.  Fiel  d  of  vuew  20  nun.  Frame  t  iiitefi/ub :  a,0;b,l4. 

10.  A  PP  disc  (light)  pressed  into  an  HMX  pellet.  The  first  frame  is  before 
impact.  Considerably  bulk  deformation  of  the  I’P  occurs  but  there  is 

no  ignition  of  the  explosive.  Field  of  view  20  mm,  Frame  time6/ps: 
a , 0 ; b , 9 1 ; e , I A 1 ; d , 2 1 4 ; e , 247 ; f , 402 . 

11.  Similar  to  figure  It)  but  wild  a  PC  disc;  except  that  impact  lias  already 

occurred  in  frame  l.  The  PC  fait;,  eat  astrophi rally .  The  first  of  several 
ignition  sites  can  bo  deleted  in  f » amc  d.  Field  of  view  20  mm.  Frame 
time/ ns  ;u,(J  ;b ,  34  ;e  ,94  ;d ,  I  2  I  1 .15  ;  f  ,  147. 
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Experimental  arrangement  for  obtaining  pretisure-t irne.  curver,.  W,  drop- 
weight;  H,  Rl,R2,lO,  hard  steel  rollers;  C,  cylindrical  guiding  sleeve 
S,  sample. 

Stress-strain  curves  for  samples  of  (a)  PP,  (b)  PC.  Catastrophic 
failure  of  the  PC  sample  occurs  at  a  strain  of  ca.  1.8.  No  discontin¬ 
uities  are  observed  in  the  PP  curve. 

Experimental  arrangement  for  friction  experiments.  D,  detector; 

S,  slits;  CD,  chopper  disc;  SD,  sapphire  disc;  PP  polymer  pin; 

BA,  balance  arm;  W,  counterweight;  DS,  drive  shaft. 

Calibration  curves  for  the  friction  apparatus  (a)  high  temperature 
region  (b)  low  temperature  region. 

Examples  of  the  traces  obtained  with  the  apparatus  (a)  is  a  soldering 
iron  at  575  k  (b)  a  small  diameter  PC  pin  and  (c)  a  large  diameter 
PC  pin.  H  is  the  signal  strength  when  there  is  no  filter  between 
source  and  detector  and  H*  the  strength  when  a.  PMMA  filter  is  inter¬ 
posed  between  source  and  detector. 

Polymer  pins  after  friction  experiments,  (a)  a  PC  pin  showing  many 
layers  which  have  been  sheared  off  and  built  up  on  the  trailing  edge 
(b)  a  PP  pin  which  has  rounded  end  shows  evidence  of  melting  (c)  a 
PMMA  pin  which  haa  a  roughened  surface  caused  by  small  port  Idas 
being  brittly  chipped  out  of  the  surface. 

Micrographs  of  polymer  pin  surfaces  after  friction  experiments. 

(a)  PC  and  shows  considerably  strain  in  tha  frictional,  direction. 

(b)  PMMA  which  has  a  powdery  appearance  with  many  small  particles 
chipped  out  of  the  surface. 
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FIG  13 


The  remaining  three  sections  of  this  report  (3,4,5)  are  on  reaction 
kinetics  and  their  analysis.  All  have  been  published  in  the  open  literature. 

For  complateneaa,  reference  is  given  to  a  paper  (Chaudhri  et  al,  1977) 
baaed  on  work  performed  by  Dr.  Chaudhri  when  he  visited  Picatinny,  and  to 
other  papers  published  by  the  group  since  1977. 
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from  Energetic  Materials,  Vol.  I,  Ed.  H.D.  Fair  and  R.F.  Walker  (Plenum, 

1977). 

J.T.  Hagan  and  M.M.  Chaudhri,  "Fracture  Surface  Energy  of  High  Explosives 
PETN  and  RDX",  J.  Mater.  Sci.  \2  (1977)  Letters  1055-58. 

M.M.  Chaudhri,  W.L.  Garrett,  0.  Sandus  and  N.  Slagg,  "The  High  Velocity 
Detonation  of  Single  Crystals  of  Alpha-Lead  Aside",  Propellants  and  Explosives, 
2,  91-93  (1977). 
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(1978).  ~ 

T.B.  Tang  and  M.M.  Chaudhri,  "Dielectric  Breakdown  of  Ionic  Azides",  Conf. 
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AltS  I  K  AC V|' 

The  paper  describes  a  graphical  computer  method  for  analysing  TO  and  DSC 
truces  which  gives  nil  three  reaction  parameters  </■',  A,  and  n)  characterising  an  mil- 
order  reaction  from  a  single  truce.  If  the  reaction  proceeds  in  multiple  stages  (£',  A 
and  n)  triplets  can  be  obtained  for  each  stage.  The  computer  programmes  ore  basically 
simple  and  use  little  computing  time  (typically  a  few  seconds).  The  advantages  of  the 
approach  over  earlier  methods  are  discussed.  As  a  test  of  the  method,  results  on  the 
dehydration  of  calcium  oxalate  monohydrate  nre  described. 

INTRODUCTION 

The  deduction  of  the  reas  on  parameters  E  (activation  energy),  A  (frequency 
factor),  and  n  (reaction  order)  from  TG.  and  DSC  traces  is  clearly  of  interest,  A 
number  of  methods  have  been  used  in  the  past  to  do  this  but  all  have  unsatisfactory 
features.  In  the  work  described  here  a  number  of  Fortran  programmes  were  written 
to  extract  all  three  reaction  parameters  of  an  mh-ordcr  reaction  from  a  single  TG,  or 
DSC  trace.  The  programmes  were  run  on  the  Cambridge  IBM  370,165  computer, 
and  the  plotter  was  used  for  graphical  output.  So  that  the  advantages  of  the  present 
approach  can  be  appreciated  earlier  methods  arc  briefly  discussed. 

HARUIR  Mf.rtlODS 

Murray  and  White1  have  shown  that  the  temperature  at  the  reaction  peak  of 
DTA  traces  cun  be  used  to  obtain  E  and  A  from  a  number  of  traces  if  a  first  order 
•reaction  (cqn.  I)  is  assumed,  where  u  is  the  fractional  residual  weight,  t  the  time,  T 
the  temperature  and  R  the  gas  constant, 


If  cqn.  (1)  is  differentiated  to  obtain  the  maximum  for  dir, 'dr  then  cqn.  (2)  can 
be  derived,  where  </>  >»  dT/dt  i>  the  heating  rate  and  Tm  the  temperature  at  the  rule 
maximum. 


A  plot  of  l/Tm  versus  in  W»/r,J)  yields  E/R  as  the  slope  and  In  (AR,'E)  as  the  intercept 
of  the  straight  line  with  the  ordinate,  The  main  disadvantage  of  this  method  is  that 
it  only  uses  a  single  point  (T„)  from  a  trace  to  produce  one  point  on  the  activation 
energy  plot.  A  number  of  experiments  urc  needed  to  obtain  a  reliable  value  for  £. 
There  is  also  the  limitation  that  the  theory  only  works  for  a  fust-order  reaction. 


Kissinger''  shnwei!  that  the  tc(.liiii(;>ii*  can  he  extended  loan  mii-oider  rent  lion 
by  deriving  cqn.  O). 


~  A  n  w",“ 1  o;p 
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He  shows  in  his  paper  that  Ihc  product  n  w|,n'1  is  independent  of  the  healing 
rale  (i/.i)  and  almost,  equal  to  unity.  Thcrv.foit,  cqn.  (.’<)  effectively  i  educes  to  vqn.  12), 
The  approach  again  makes  poor  use  of  the  information  contained  in  n  DTA  truce. 
Ottawa's  analysis''  of  TC  traces  begins  with  a  general  form  of  reaction  equation, 
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df 
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(4) 


This  is  intcgiatcd  by  Ozawa  wlio  shows  that  the  solution  to  the  right-hand  side  of  the 
equation  is  a  polynomial  in  the  variable  fi/RT,  liquation  (5)  is  an  approximate 
formula  which  relates  the  heating  rate  <f>  with  /■/. 


log  <j>x  +  0,4567  -£.~  «  log  4>,  +  0.4567  ~  (5) 

where  7\  and  T2  coricspond  to  the  same  fractional  decomposition  («.')  of  the  sample. 
Since  a  plot  of  log  4>  versus  I /Tcan  be  obtained  for  ouch  value  of  w,  Ozawa  aggregates 
all  these  results  in  a  master-curve  with  a  gain  in  accuracy  of  the  final  result,  liquation 
(5)  only  assumes  thut  the  function  f(iv)  is  constant  for  n  given  value  of  w. 

This  method  also  needs  a  number  of  experiments  ut  different  heating  rates  to 
generate  enough  data  for  a  single  activation  energy  plot.  Although  the  method  was 
developed  for  thermogravimeuie  data,  it  can  be  applied  to  DSC  data  by  integrating 
the  peaks  to  give  values  of  »\  A  theoretical  analysis  of  such  a  process  is  given  by  Reed 
cl  al.4. 

Borelwrdt  and  Daniels*  describe  n  method  which  allows  E,  /<  and  n  to  be 
determined  from  the  shape  of  a  DTA  curve.  Their  approach  was  applied  to  reactions 
in  the  liquid  phase.  'I  lie  type  of  reaction  to  which  this  method  is  applicable  must  have 
n  tingle  rule  constant  and  the  activation  energy  most  not  vary  with  temperature. 
By  considering  the  heat  transfer  equation  of  the  calorimeter  cqn.  f6)  was  obtained 
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where  A'  is  the  number  of  moles  present,  ,VU  the  initial  number,  A  the  heat  transfer 
coefficient,  .<  the  area  under  the  I)  I  A  curve,  C,.  the  heal  capacities  of  the  Luo  liquids 
and  ,i  7  the  lempeiature  diilercnce  between  ihe  two  cells.  I  mediation  of  cqn.  (6)  wuh 
respect  to  time  gives 
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where  a  is  the  area  under  the  DfA  curve  up  to  the  present  temperature. 
The  cxpic.vsiori  tor  the  rale  constant  k  of  a  reaction  of  order  n  is 


k  rr  .  ..  y*m~ 11  /d< 

N* 

where  V  is  the  volmm-.  Substituting  cqns.  (6)  an  I  (7)  for  d.V/dr  and  A  gives 


(8) 
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In  order  to  olilnin  reaction  p  immoier .,  Uou'li. licit  ami  (,>aiiul-‘  assumed  n  wine 
for  ii.  This  allowed  k  to  be  inhul.iled.  A  plot  of  In  k  xeiMiS  I//  (;u  fixation  Cilia  ry 
plot)  yields  a  Mraiyht  line  if.  and  only  if.  the  cm  reel  \;dne  of  i,  liar,  been  assumed, 
tints  valttcs  of  A  and  n  can  be  obtained.  The  main  dillicviliy  of  this  method  is  the 
number  of  plots  that  lttivu  to  be  made  for  cadi  value  of  u  until  the  points  lie  on  a 
straight  line. 

A  number  of  computer  programmes  are  presented  by  Itenin  ct  al.6  and  used  to 
solve  three  problems,  namely  (i)  the  production  of  thermograms  given  the  reaction 
parameters  Ji,  A,  n  und  given  that  the  reaction  is  of  the  /ith-order  type;  (ii)  the  calcula¬ 
tion  of  the  (/;,  A,  n)  triplet  from  an  experimental  thermogram  this  is  the  inverse  to 
problem  (i);  und  (iii)  the  determination  of  the  reaction  parameters  and  the  complex 
reaction  mechanism  from  a  given  experimental  thermogram.  The  first  problem  is  a 
st ru I ght forward  function  plot  of  the  solution  to  the  wilt-order  reaction  equation. 
Tito  second  problem  is  more  difficult,  and  presents  itself  to  anyone  who  wants  to 
obtain  reaction  parameters  front  experimental  data.  Benin  ct  al.*  treat  it  as  a  problem 
in  curve  At  ting.  They  minimise  the  discrepancy  between  the  experimental  data  and  the 
calculated  function  by  taking  as  the  discrepancy  criterion  the  maximum  deviation  of 
the  function  from  the  data.  The  solution  of  the  third  problem  Is  an  attempt  to  reveal 
the  mechanism  of  a  chemical  reaction.  The  programmes  become  quite  involved  as 
they  work  ottt  activation  energies  for  each  section  of  the  thermogram  by  using  a  plot 
or  reaction  rate  versus  1/7  for  a  constant  amount  of  con, version,  a,  but  different 
heating  rates. 

The  most  general  form  of  the  (It  function  assumes  a  decomposition  process 
which  proceeds  In  a  numboi  of  individual  stages.  Kadi  stage  consists  of  an  elementary 
nth-order  reaction,  i.c,,cqn.  (I )  with  A,  /-.'and  ir  replaced  by  Ar  Et  and  u'i,  respectively, 
whore  *7”  denotes  the  vuriable  belonging  to  the  /tit  stage  of  the  reaction,  The  program¬ 
me  which  opllmisestho  At  of  this  function  to  the  data  lakes  up  to  4  h  of  computing  time. 
This  make*  the  method  rather  expensive.  Another  weak  point  Is  the  discrepancy 
criteria  between  the  data  and  the  lit,  which  is  taken  as  the  maximum  difference,  It  is 
well  known  that  the  data  points  towards  the  end  of  a  decomposition  process  become 
less  reliable  due  to  impurities  and  zero  errors.  Hut  these  points  have  the  same  weight 
in  the  analysis  as  all  the  others. 

Schempf  ct  ah’  describe  a  programme  which  fits  a  least  square  polynomial  to 
the  weight-loss  curve  The  quantity  du/df  is  calculated  by  dilTcreniiatiun  of  the 
polynomial.  A  straight  line  is  then  titled  to  it  plot  of  log  k  versus  1/7'.  and  the  least 
square  ciitcrion  is  used  again.  The  reaction  order  is  assumed  to  be  one,  and  nd 
attempt  is  made  to  eliminate  the  deviations  from  an  ideal  reaction  ut  the  beginning 
and  end  of  the  weight-loss  curve.  Out  work  shows  that  this  is  very  critical  and  enn  be 
obtained  from  an  Arrhenius  plot  performed  by  the  computer. 

NliW  METHOD 

The  method  developed  was  based  on  the  ideas  of  llorchardt  and  Daniels',  Our 
data  were  obtained  from  both  a  thcrmoh.ikmcc  (.Stanton  RcJergt't)  and  a  differential 
thermal  calorimeter  (Pci ki n-  [  liner  DSC'-."’  I 

The  analysis  assumes  a  reaction  which  can  r*c  described  by 


"  -  At  exp  (  -  -  ®J.  )  f(ii-) 


Since  the  heating  rate  h  (deg  min  '1)  is  cvnxlmt,  the  time  variable  can  be 
eliminated  ft  out  cqn.  (10)  by  punim'.  d  r  dr  1  h/OU  (d,\v.  see  ' 


Taking  natural  lojwi  iihni.v  on  both  sidv ■■,  fives 


(U) 


Note  that  (••■  dt»/d/)  is  ti  positive  quantity  since  the  Frnctionnl  residual  weight  decreases 
with  increasing  temperature. 

A  plot  of  {ln(-  die, Air)  -  In  f,(»r)}  versus  (I (T)  only  yields  n  struigln  line  if 
fiO»')  has  been  chosen  correctly.  In  the  case  of /,(tt‘)  «  n"  the  reaction  order  it  is 
varied  tin'll  n  straight  line  Is  found,  The  activation  energy  E  con  be  obtained  from  the 
slope  of  this  plot,  the  frequency  factor  A  from  the  intercept  on  the  ordinate  and  the 
reaction  order  n  is  determined  by  the  straight  line  criterion. 

Tills  analysis  can  also  be  applied  to  the  DSC  data  since  the  basic  eqn.  (12)  only 
differs  from  (10)  by  constants. 


~j$‘  "  ±  AHty,  “  "-,)  A  exp  ^  —  f(w)  (12) 

The  +  (--)  its  the  equation  corresponds  to  an  exothermic  (endothermic) 
reaction  where,  Q  is  the  hunt  given  off  by  the  sample,  Alt  is  the  enthalpy  change  of  the 

entire  reaction,  »>0  is  (ho  initial  weight  and  tr,  is  tiie  rest  weight,  and  the  other  para¬ 
meters  arc  us  defined  previously. 

Taking  logarithms  gives 
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This  time  dg/dl  at  temperature  T  is  the  measured  quantity  and  »•  needs  to  bo 

calculated.  This  is  easily  done  since,  if  the  reaction  results  in  n  weight  loss 

» 


r* 


The  computer  numerically  integrates  d Q}\\t  to  yield  the  weight  at  time  l.  A  plot 
of  {In  <\Qldt  -  la  f(a  )}  versus  1/7  produces  l',  A  and  u  (if  f(w)  «>  u^Jin  Ihc  same  way 
as  for  the  weight- loss  curves  above.  The  slope  is  again  -£//?  and  the  intercept  is 
In  d//(w0  -  w ,)A. 

In  practice  the  ditto  from  the  TO  or  DSC  is  recorded  on  8  track  ASCI  paper 
tape.  This  Information  is  then  rend  into  dive  files  c»T  the  computer.  These  files  can  be 
used  by  the  programmes  which  calculate  the  activation  energy,  etc.  After  these 
programmes  have  i  ua,  tltc  plotter  output  is  viewed  on  o  television  screen.  Only  if  the 
plot  is  satisfactory  is  u  copy  obtained  from  the  plotter  for  the  final  analysis.  The 
dimensions  correspond  to  the  original  output  of  the  chart  recorder  so  thnt  a  direct 
comparison  is  possible.  Different  pen  colours  nre  used  to  differentiate  between  the 
nuMlntn  (black)  and  the  theoretical  fit  to  the  data  (green),  these  colours  are  re¬ 
presented  in  Fig.  1  ns  a  broken  line,  (black)  nod  n  full  line  (green). 
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Fifi.  1.  Weight  low  curve  or  Ca  C;0<H?0  analysis.  The  dashed  line  is  ihc  original  data  and  the  full 
line  Ihc  replot  using  the  computed  ( F.,A,u )  triplet  for  the  three  straight  line  sections  which  arc 
separated  by  vertical  dashed  lines. 


The  differential  programme* 

There  are  three  different  ways  in  which  differentials  to  the  weight-loss  curve 
were  obtained.  The  ru  st  uses  differences  and  produces  a  noisy  differential  curve,  which 
is,  however,  nearest  to  the  original  data.  The  other  two  use  a  polynomial  fit,  and  a 
spline  fit,  respectively.  The  polynomial  fit,  usually  or  25th  order,  is  generally  good  but 
produces  artificially  large  undulations  over  straight  regions  of  the  curve.  The  best 
results  were  obtained  by  a  cubic  spline  fit.  This  programme  (its  a  third  order  curve  to 
sections  of  the  data,  and  joins  them  in  such  a  way  that  the  function  and  the  first  two 
derivatives  arc  continuous.  The  programme  uses  a  standard  subroutine  from  the 
Harwell  Fortran  library,  which  chooses  its  own  knots  (points  where  two  splines  meet) 
according  to  the  behaviour  of  tiic  original  curve.  This  means  that  over  a  straight  part 
of  the  curve  a  long  section  can  be  fitted  by  a  single  cubic  spline,  whereas  at  places  of 
greater  variation,  smaller  sections  are  chosen  to  be  approximated.  The  differential 
of  the  curve  is  then  obtained  by  differentiating  the  third  order  spline. The  fit  of  the 
spline  and  the  differential  is  very  accurate.  The  amount  of  data  smoothing  achieved 
by  this  fit  suits  the  purpose  of  the  following  programme  (next  section)  best. 

The  programme  to  find  E,  A  and  it 

ll  is  assumed  in  this  analysis  that  the  reaction  that  leads  to  the  weight-loss  curve 
is  at  least  piece-wise  of  the  nth  order,  i.c.  follows  eqn.  (10)  with  f;(tv)  =  This 
programme  plots  (In  (-  dtv/dr)  -  «  In  tv}  versus  \IT  for  different  values  of  n.  The 
values  of  du/d/  and  of  tv  arc  obtained  from  a  file  that  was  created  by  the  previous 
programme  and  contains  the  smoothed  weight-loss  data  and  the  differential  in  the 
form  of  the  spline  parameters.  The  value  of  n  is  taken  as  zero  for  the  first  plot  and 
then  increased  by  0.2  let)  limes  up  to  n  «•  2.  These  eleven  lines  have  different  curvatures. 
The  straightest  of  them  can  easily  be  selected  by  eye  and  this  n  value  is  then  taken  as 
the  reaction  order.  The  slope  yields  F.  and  the  intercept  A.  The  resulting  triplet  or 
numbers  (£,  A,  a)  are  necessary  and  sufficient  to  characterise  an  nth-order  reaction. 
The  type  of  reaction,  i.c.  f(iv)  tv"  can  be  changed  easily  so  that  the  programme  is 
applicable  to  reactions  of  any  type,  i.c.  f(tv)  can  have  any  form  (e.g.  random  chain 
scission). 

*  This  mu!  the  following  pn  >jt  a  mines  can  be  obtained  from  the  authors. 


32. 


The  programme  to  produce  a  theoretical  Jit  from  the  (F.t  A,  n)  triplet 

This  pronrammc  allows  the  production  of  a  weight-loss  curve  with  up  to  four 
different  sections  each  of  which  has  its  own  E.  A  and  «  factor.  It  was  implemented  on 
a  Hewlitt  Packard  Desk  top  calculator  and  plotter  of  the  type  HP  9S25.  The  theoretical 
plots  produced  in  such  a  way  fit  the  original  weight-loss  curve  very  well.  They  also 
provide  a  useful  check  that  all  the  calculations  necessary  to  produce  this  triplet  were 
correct,  and  produce  a  good  fit. 


V 


Fig.  a.  Activation  energy  plot  for  various  values  of  reaction  order  n.  The  «  value  which  gives  live 
best  straight  line  segments  is  the  one  appropriate  to  the  reaction  (here  0.6).  There  arc  three  steps  to 
the  reaction  but  the  fust  two  correspond  to  only  ca.  2%  of  the  weight-loss. 


TEST  OF  METHOD 

The  dehydration  reaction  of  calcium  oxalate  monohydrate 
Ca  CjO^O  -►  Ca  C204  -f  HaO 

is  well  understood  and  yalucs  for  the  reaction  parameters  can  be  obtained  from  the 
literature8,  ’.  It  therefore  provides  a  good  test  ease  for  the  new  analysis  method. 
Figure  1  shows  the  weight-loss  curve  of  a  sample  of  Ca  Cj04H,0.  The  dashed  curve 
represents  the  original  data  and  the  full  curve  the  theoretical  fit  based  on  the  computed 
( E ,  A,  n)  triplet.  Figure  2  shows  the  activation  energy  plot.  The  number  at  the  end  of 
each  line  gives  the  reaction  order  that  is  assumed  for  the  calculation  of  the  corre¬ 
sponding  line.  The  straight  lines  which  approximate  the  curves  best  arc  used  to  obtain 
(£,  A,  n)  triplets  for  the  three  sections.  Since  these  (E.  .-I.  n)  triplets  compleicly 
specify  the  reaction  the  original  weight-loss  curve  can  be  modelled  (full  curve  Frig.  I ). 
The  results  obtained  with  this  method,  and  literature  values  are  summarized  in  Table 
1.  AH  our  experiments  were  performed  at  a  heating  rate  of  10  deg  min"1.  The  results 
obtained  in  this  way  have  an  error  of  ±  5  kJ  mol" 1  which  is  due  to  the  variance  of  the 
samples  themselves.  The  method  itself  lias  much  greater  accuracy.  It  was  found  that 
the  best  straight  line  could  be  filled  to  a  curve  corresponding  to  a  reaction  order  of  0.6 


TABLE  1 


E  hi  A  Kef. 

(kJ  mole'*) 


87 

19.8 

Tills  work 

14 

17.5 

This  work 

as 

21.0 

X 

n 

9 
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which  indicates  a  surface  reaction  (theoretically  2/3).  It  can  also  he  seen  on  Fig.  2 
that  there  arc  three  distinct  stages  of  the  reaction.  The  first  two.  however,  correspond 
to  only  2%  of  the  weight  and  arc  attributed  to  volatile  substances  on  the  surface. 
This  shows  the  sensitivity  of  the  method  used  especially  considering  that  the  vertical 
axis  has  a  logarithmic  scale.  The  real  dehydration  reaction  is  represented  by  the  last 
straight  line  segment. 

The  method  outlined  in  this  paper  has  recently  been  successfully  applied10  to 
decomposition  studies  of  P.E.T.N.,  P.E.T.N.  with  fillers,  high  density  polyethylene 
(HOPE)  and  polytctrafluoroclhylcnc  (PTFE). 
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In  most  solid  state  reaction*  the  reaction  velocity  can  be  described  as  a  product  of 
two  functions  KlT\  and  /<1  —  »)  where  7  is  the  temperature  and  *  the  degree  of  con¬ 
version  of  the  solid  reactant.  The  ph>sieal  interpretation  of  these  functions  is  discussed, 
and  a  systematic  metnod  is  described  by  which  At  -  a)  of  a  reaction  is  identified  from 
its  kinetic  data.  K{T)  and  the  reaction  mechanism  are  then  determined.  This  method 
has  been  successfully  applied  to  analyse  the  kinetics  of  the  thermal  decomposition  of 
silver  azide. 

In  a  solkl-statc  reaction,  the  reaction  velocity  is  given  by  —  d/dr  (i  -  a)  ■>  «, 
where  a  *  a (/)  is  the  fraction  of  the  solid  reactant  which  has  reacted  by  time  f. 
Its  kinetics  can  be  solved  by  determining  a  ns  a  function  of  temperature  and  the 
global  amount  of  reactant  left.  This  phenomenological  knowledge  is  n  necessary, 
though  not  sufficient,  condition  for  elucidating  the  reaction  mechanism.  It  is 
necessary  to  formulate  the  reaction  velocity  in  terms  of  the  global  variable  a, 
because  there  is  a  continuous  collapse  of  structure  in  the  reactant.  Furthermore, 
the  local  concentration  of  the  reactant  varies  throughout  the  reaction  volume  and 
cannot  be  used  as  a  state  variable.  In  fact,  unlike  the  ease  of  a  homogeneous  re¬ 
action  in  the  liquid  or  gaseous  phase,  there  is  no  real  ’reaction  order’  with  respect 
to  any  reactant  in  a  reaction  involving  condensed  matter  whose  mechanism  is 
usually  of  the  heterogeneous  type. 

If  the  reaction  proceeds  isothcrmally,  it  is  observed  empirically  that  (a.  i) 
curves  corresponding  to  different  temperatures  T  arc  isomorphic  to  one  another, 
at  least  within  a  range  of  T,  i.e.,  by  a  linear  scale  chance  in  t.  different  curves  can  be 
superimposed  [l].  It  follows  that  i  is  a  separable  function: 

(a}iSo,b«n-,..  =  TO/? I  -  *)•  (I) 

Hcrc/(1  -  a)  may  change  in  different  ranges  of  T  or  x  For  every  /{I  -  a),  there 
corresponds  a  single  A '(7).  It  should  be  noted  that  experimental  data  may  be  ade¬ 
quately  analysed  by  (1)  only  if  it  has  been  ensured  that  the  temperature  distribution 
in  the  sample  is  sufficiently  uniform  and  constant.  Furthermore,  the  theoretical 
significance  of  A(T)and  /(I  -  ar)determined  from  the  data  should  always  be  exam¬ 
ined  with|  rc2ard  tof  the  class  of  mechanisms  they  indicate.  A  question  of  con- 
islency  arises  in  this  respect.  In  the  literature,  coherent  and  integrated  accounts 
of  the  physical  (as  contrasted  with  the  Formal)  meanings  of  K(T)  and  /(I  —  x) 
are  not  easily  found.  A  discussion  on  their  interpretations  therefore  forms  the 
first  part  of  this  paper.  In  the  second  part  we  propose  an  efficient  method  of  de¬ 
termining.  with  confidence,  both/(l  •—  a)  and  A"(7*)  from  0  and  z(/)data.  The  inter¬ 
pretation  of /(I  -  x)  is  the  essential  link  in  this  method  between  the  experimental 
data  and  the  functions. 


Physical  interpretations 

The  function  f( /  —u) 

Solid-state  reactions  arc  com  -rex  processes  which  proceed  in  several  stages. 
These  can  be  the  delocalization  >r  transfer  of  electrons  in  chemical  bonds  (in  the 
case  of  non-metals),  (he  diluisio  i  of  aioras.  free  radicals,  or  ions.  Ihc  desorption 
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of  product  molecules  when  they  are  in  the  gaseous  phuse.  the  heat  transfer  to  the 
reuction  zone  in  the  case  of  endothermic  reactions,  and  the  formation  of  a  ntw  solid 
structure  (crystalline  or  amorphous)  if  one  of  the  products  is  in  the  solid  phase. 
The  last  step  may  often  be  further  differentiated  into  nucleating  growth  of  nuclei 
(at  velocities  which  depend  on  sizes  of  nuclei)  [i )  and  sometimes  the  collapse  of 
the  lattice  from  a  transitory  one  to  the  equilibrium  structure  [2].  Irrespective 
of  the  details  of  reaction  mechanisms,  however,  under  a  given  set  of  circumstances 
(T.  a,  sample  history  etc.)  one  of  the  stuges  will  be  the  slowest.  H  then  acts  as  the 
rate'  limiting  step  cf  the  reuction.  and  it  will  determine  the  kinetics  i.e..  the  rate 
luw  in  (I). 

Further,  a  solid-state  reaction  has.  in  contrast  to  a  heniogertous  reaction  whose 
progress  is  independent  of  spatial  coordinates,  an  additional  controlling  factor, 
numely  lopochemistry.  This  refers  to  (he  geometrical  shape  of  the  .solid  reactoot 
and,  in  different  cases,  to  its  free  surface  area,  its  defect  structure,  the  thickness 
of  the  product  layer  if  solid,  or  to  the  product- reactant  boundary,  etc, 

The  function /U  -  *)  reflects  the  nature  of  the  rate-limiting  step  and  the  topo- 
chemistry  of  the  reuction.  It  muy  accordingly  depend  on  certain  sample  conditions, 
such  as  whether  the  sample  has  been  pre-irradiuted  or  bleached,  and  whether 
the  sample  Is  in  the  form  of  a  powder,  or  a  large  single-crystal  of  a  different  shape 
from  the  crystallites.  It  will  vary  in  several  ranges  of  T  If  in  each  of  them  a  difterent 
elementary  step  becomes  rate-limiting,  as  occurs  in  the  decomposition  of  potassium 
azide  [3],  At  u  given  temperature,  it  may  also  change  In  different  ranges  of  a. 
due  to  the  switching  of  the  rate-limiting  step  or  topochemical  changes.  This  hap¬ 
pens,  for  instance,  in  the  oxidation  of  zirconium  [4]  and  in  most  decomposition 
processes  [1].  An  extreme  case  is  the  decomposition  of  ammonium  perchlorate 
which,  in  the  two  temperature  regimes  below  and  above  620AT,  has  entirely  dif¬ 
ferent  reaction  mechanisms  and  in  fact  yields  different  reaction  products  [5]. 
In  such  cases  there  may  be  competing  paths  for  the  chemical  reaction  or  it  may 
in  fact  be  followed  by  another  chemical  reaction  whose  'onset  temperature’  is 
higher.  In  all  possibilities,  however,./?!  -  a)  should  be  the  same  for  u  given  a 

independent  of  /’  (within  a  range)  if  it  Is  to  have  more  than  only  an  empirical 
sign  iflea  nee. 

In  Table  1  we  have  collected  together  the  more  common  forms  of  ./?  I  -  a) 
which  have  been  used  in  the  literature,  and  the  corresponding  integrated  forms 

[  ■ 

\  d /  -  a)  -  |  da//?l  -  a)  s  F(v).  Note  that  F(a)  »  A  U  -  /„)  if  the  runge 

o  o 

of  a  for  which  it  becomes  applicable  starts  at  a  -  a (/„).  Also,  for  simplicity  here¬ 
after  we  write  K  for  K(T), 

In  many  reactions,  such  as  most  dtccrrf  esitiers  and  dehydrutions,  the  rate- 
limiting  step  takes  place  at  the  inter!) cc  betwee  n  different  phases  as  in  sublimation. 
The  speed  at  which  the  interface  rr.uvu  d.w  .suclant  is  (ut  n  given  temperature) 
then  either  a  constant,  or  a  unique  function  of  the  interfneial  area.  This  area 
therefore,  from  the  kinetic  point  of  view,  plays  the  same  role  as  that  of  con¬ 
centration  in  homogeneous  reactions.  If  the  speed  is  constant,  then  the  theoretical 
significance  of  /? I  -  a)  is  blear:  it  gives  the  areu  expressed  us  a  fraction  of  the 
original  area  at  a  ■  0.  This  is  the  case  of  a  reaction  controlled  by  the  movement 
of  a  coherent  phase-boundary  and  listed  us  t\  G  and  H  in  Table  I.  In  this  situation, 
the  explicit  form  of  ./?  I  -  a)  depends  on  the  geometry  of  the  reacting  system, 
though  generally  it  is  a  decreasing  function  of  x  or  at  most  constant. 

If  the  reaction  consists  of  the  formation  of  compact  nuclei  of  a  solid  product 
at  localized  places  in  the  rcoctuni  followed  by  their  relatively  rapid  growth,  then, 
to  express  the  total  interfaciul  areu,./?l  -  a)  is  derived  from  the  luws  of  nucleatlon 
and  growth.  This  is  the  situation  when  the  reaction  is  uutocutulytic  (13):  reactant 
molecules  ut  u  reactant-product  interface  react  in  preference  to  those  at  a  reactunt- 
’vacuum'  surface.  The  preference  is  due  to  the  existence  of  microstrains  in  the 
reactant  at  the  Interface,  or  due  to  the  electrochemical  potential  of  the  product 
phase  when  the  rate-limiting  step  is  a  redox  process.  The  various  possible  forms  of 
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/'( I  -  for  ;in  autocatulylic  reaction  arc  lifted  in  Table  I .  A  to  F.  Noie  that  they 
have  the  general  form  /'(I  -  arl  -  I  -  vi11.  In  A.  F  and  £.  </  is  zero  and  i 
increases  monotonicallv  with  x.  Such  a  situation  is  most  unlikely  to  iust  up  to 
y  *  1.  These  types  of. /(I  -  *)  therefore  may  apply  only  to  the  acceleratory  purt, 
if  it  is  present  in  the  a-tinie  curve  and  which  will  usually  be  followed  by  a  decay 
pun.  In  C  and  D,  q  is  non-zero  and  these  types  give  sigmoid-shaped  curves.  The 
inflexion  point  occurs  at  a'  =>  p!(p  +  q),  as  can  be  seen  at  once  from  the  condition 
i<  =>  0.  It  is  thus  kinetically  feasible  for  them  to  lit  the  complete  experimental  curve. 
There  are.  however,  physical  grounds  to  consider  that  even  they  should  be  used 
to  analyze  only  the  accelerutory  period  [I], 

In  other  reactions  the  rate-limiting  step  is  not  eonflned  to.  or  does  not  only 
occur  nt.  the  rcuctant  surface.  For  instance,  in  the  imimolecular-tlecay  type  or 
reaction,  all  molecules  whether  on  the  surface  or  in  the  bulk  huve  an  equal  proba¬ 
bility  per  unit  time  of  reacting.  This  is  the  case  when  the  change  from  the  re¬ 
actant  to  the  solid  product  phase  involves  little  re-arrangement  of  the  reactant 
atoms.  The  reaction  has  homogeneous  mechanism  and  thus  a  true  reaction 
order  of  one,  i.e. /II  -  x)  «■  I  ~  *.  Many  decomposition  reactions  tend  to  this 

Table  I 

The  common  types  of  .solid-state  reaction 


Ruction 


Autocatalytlc  Power-law  nucleation  and  growth  at  constant  speed:  A 

Linear  branching  chain  of  nuclei,  no  overlap  during  growth:  B 
Branching-chain  nucleation,  interference  during  growth:  C 
Random  nucleation,  growth  accompanied  by  ingestion  of  nuclei :  D 
Instantaneous  nucleation.  size-dependent  growth:  E 


Phase-boundary  controlled  decay  1  -dimensional :  f 


Unlmolecular  decay:  / 

Diffusion-controlled 

yxt  -  «)  -  an 

2- dlmenslonal:  G 

3- dlmensional :  H 

1  -dimensional :  J 

2- dlmensional:  K 

3- dlmensional:  L 
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•*-  i 

i)  Log  {—In  (I  —  *)[•  as  Constant  +  in  log  r.  0.15  <  *  <  0.5 

il)  Plot  of  L.H.S.  in  i)  against  log  r  distinctively  concave  upwards;  f  „  and  r„  arc  constants 
while  C(D  is  function  of  temperature 

ill)  h  -  0.774,  0.700.  0.664,  0.642 .  0.356  for  rt  *  2.  3.  4.  5.  ...  75 

Is)  Alternative  derivations;  we  have  obtained  the  appro  si  mate  forms  of  f(  1  -  i)  by  es- 
pandlng  into  series  (1  -  a)-'1,3  and  (1  -  a)'1'3  to  second  order  in  /:  resulting  error  st  **/ 6 
l<  10 "i  for  r  <  0,8) 


limit  at  high  a  values.  Another  example  is  when  the  rate-limiting  step  is  the 
migration  of  product  ions  along  the  dislocation  network  to  form  additional 
growth  nuclei  at  dislocation  nodes  f  17].  This  leads  to/T  1  —  ot)  =  ot,  the  same  form 
as  for  branching  nuclei  [7J.  The  slow  process  of  ‘ageing’  in  some  explosives  when 
they  are  stored  at  room  temperature  may  be  by  such  a  nucleus-chain  mechanism. 
A  third  category  comistc  of  reactions  controlled  by  the  diffusion  of  reactants 
across  a  product  lajJer  which  is  solid.  The  diffusion  may  proceed  uniformly  through 
the  bulk  of  the  layer  and  is  thus  sLucture-inscnsitive,  or  preferentially  along  its 
gross  lattice  imperfections  arising  from  product-reactant  mismatch,  in  the  case 
of  uniform  diffusion,  the  speed  at  which  the  produci-reactant  interface  moves  is 
u  function  only  of  the  product  thickness  (and  temperature),  and  the  appropriate 
forms  of f{\  -  «)  are  included  in  Table  1  as  J,  K  and  L,  The  oxidation  of  metals 
often  follows  diffusion-controlled  kinetics;  in  sheet  form  these  tarnish  according 
to  the  parabolic  law  ct*oc /.  Exceptions  uret  hose  metals  in  Groups  la  and  Ila  of 
the  Periodic  Table.  Excluding  beryllium,  they  all  form  oxide  layers  which  are 
porous,  so  that  the  atoms  of  the  metal  do  not  have  to  diffuse  through  a  coherent 
layer  before  coming  into  contact  with  oxygen. 


The  Junction  K{T) 

It  is  almost  always  the  case  that  the  temperature-dependent  part  of  (I)  can  be 
represented  successfully  by :  - 

K(T)  m  K*  e\p(-EjkT)  (2) 

in  which  k  is  Boltzmann's  constant,  t’.nd  the  mncro-kir.etic  constants  £  and  K-,. 
do  not  depend  on  T (within  the  range),  though  usually  they  take  on  different  values 
when 7(1  -  «)  chunges, 

If  it  is  established  that  the  reaction  is  rate-limited  by  a  diffusion  or  migration 
process,  the  interpretation  of  K(T )  is  complicated,  but  obviously  it  is  proportional 
to  the  corresponding  transport’coefficient,  which  in  general  is  an  exponential  ftinc- 
tion  of  T.  An  over-simplified  theory  for  the  situation  of  uniform  one-dimensional 
diffusion  gives  K(T)  -  (S/V„)aD(T),  where  S  Is  the  interfacial  area,  Vn  the  initial 
volume  of  the  reactant,  and  D(T)  the  diffusion  coefficient  (cf.  [19]). 

For  a  single-solid-reactant  reaction  which  is  controlled  by  a  surface  process, 
on  the  other  hand,  the  simple  theory  of  Shannon  [18]  is  often  successful.  This 
theory  is  a  generalization  of  the  Polanyi-  Wigner  equation.  Assuming  the  existence 
of  some  activated  complex,  which  as  a  transition  state  can  be  treated  in  thermo¬ 
dynamic  equilibrium  with  the  reaotant,  he  related  the  pre-exponential  factor  K^_ 
to  the  rotational  and  other  internal  degrees  of  freedom  of  a  reactant  molecule 
in  addition  to  the  vibrational  ones.  Following  Shannon  we  can  set:  - 

-  (kT!h)  exp  (AS*fk)  6S0IV„  (3) 

and  exp  (-E/kT)  -  exp  (- AH'jkT).  (4) 


Here  the  mean-frequency  factor  kTIh  containing  Planck’s  constant  is  usually  in 
the  region  of  IOias-I(see  below),  dS*  and  d  //*  are  respectively  the  entropy  and  the 
enthalpy  of  formation  of  the  transition  complex,  6  is  the  thickness  of  one  mono- 
layer  and  y„  the  initial  volume  of  the  reactant,  and  Sn((I  -  a)  gives  the  free  sur¬ 
face  or  the  product-reactant  interface  area  when  the  degree  of  conversion  is  a. 
(Strictly  speaking,  it  has  been  assumed  that  the  reaction  proceeds  isobarically). 

Note  that  in  this  interpretation  the  empirical  quantity  Kx  contains  the  surface- 
-to-volume  ratio  and  so  depends  on  the  sample  geometry.  Also,  it  is  apparently 
proportional  to  T.  (In  gas  reactions,  the  collision  theory  gives  K,  xT*.)  In  our 
opinion,  however,  if  the  vibrational  modes  are  being  considered  then  only  at  low 
temperatures  will  the  peak  distribution  of  phonon  frequencies  lie  at  kVh.  For 
most  substances  (with  the  exceptions  of  Be,  C’r  and  diamond)  the  Debye  tempera¬ 
ture  ‘s  less  than  500K,  so  that  the  frequency  factor  should  stay  as  kO{1h  £  IOl:1s" 1 
for  all  likely  experimental  temperatures. 


The  factor  expl/1.V:  A  )  may  altcrnutively  he  wriuen  in  terms  of  partition  func¬ 
tions  as  Ql'Q.  which  cun  he  determined  from  spectroscopic  data  1 20].  In  most 
cases  cannot  be  larger  than  the  reuctant  entropy  of  melting,  and  exp (ASxjk) 
comes  out  normally  Detween  unity  und  I01.  Occasionally  explJS1'*)  is  found  to 
be  less  than  unity,  as  is  the  steric  factor  in  gas  reactions.  Such  a  negative  value  of 
d-S*  means  that  the  activated  complex  is  more  ordered  than  the  reactant  (e.g. 
[21]).  Experiments  on  some  decomposition  reactions  have  given  K  which  are  ab¬ 
normally  high  in  comparison  to  the  theoretical  values  of  (3).  Hypotheses  put  for¬ 
ward  to  explain  such  discrepancies  include  co-operative  activation  (22).  proton- 
deloculizution  (23).  and  a  mobile  layer  of  molecules  on  the  reuctant  surface  [18], 


Kinetic  analysis 

Current  practice 

In  determining  the  kinetics  one  wishes  to  find  K, .  E.  und/] I  -  st)  or  equiv¬ 
alently  F(a)  so  thut  the  reaction  velocity  cun  be  predicted  at  any  given  T  and  a. 
This  is  commonly  done  by  analyzing  a  set  of  a(M  or  equivalently  S(r)  values  ob¬ 
tained  by  monitoring  a  number  of  samples  reacting  isothermally  at  a  number  of 
temperatures.  The  consistency  of  the  Kk  and  E  values  with  /(I  -  a)  should  as  far 
as  possible  be  assessed,  and  correlated  with,  for  instance,  microscopy. 

A  quick  method  of  calculating  E  was  used  by  Haynes  and  Young  [24].  Consider 
a  set  of  (a,  t)  curves  which  have  been  found  to  be  isomorphic.  For  any  two  curves 
(a,,  /,)  und  (aj,  t./,)  corresponding  to  temperatures  T,  und  T.,  respectively,  one  can 
write 

F( a,)  -  exp  ( -£/AT|>  ^ 

F(a.;)  -  exp  ( - E/kTt ) . 

By  choosing  points  corresponding  to  the  same  a  on  the  two  curves  so  that  /'(«,)  ■ 
m  Fiit-i),  E  con  be  evaluated  by  plotting  In  r  vs.  I,T.  On  the  other  hand,  to  deter¬ 
mine  /fa)  often  a  trial-und-error  method  is  resorted  to.  Conflicting  forms  o'f  the 
function  have  sometimes  been  asserted  by  several  uuthors  for  the  same  material, 
like  NH,CI04  (sec  [25])  und  KM„0,  (see  [26]). 

A  conventional  way  of  superimposing  isothermul  curves  is  to  convert  them  into 
‘reduced-time  plots'  by  individually  scaling  their  t-axis  with  the  fuctor  1/t„  where 
t|  is  the  time  when  a  ■=  0,5  on  the  i'-th  curve.  In  this  way  A'lT,)  is  absorbed  into 
each  scale  factor  and  all  (a,  /)  «=  (0,5.  I)  points  coalesce,  while  other  *(f)  points 
may  be  plotted  out  tr  ir  the  curves  ure  indeed  isomorphic.  Shurp  et  ul.  [27] 
tabulate  the  theoretical  ues  of  t  against  i'x  for  some  of  the  F(»)  shown  in  Table  I. 
They  propose  thut  by  comparing  experimental  data  with  such  master  values  the 
correct  F(a)  can  be  identified. 

The  above  method  may.  however,  result  in  ambiguity  due  to  u  number  of  as¬ 
pects.  Experimental  data  contain  random  errors,  but  no  simple  statistical  analysis 
can  be  applied  to  the  identification  criterion  it  employs  becuuse  no  straight-line 
graphs  are  involved.  Additionally,  a  gencrul  problem  for  all  isothermal  experi¬ 
ments  is  the  zero-time  uncertainty.  The  finite  time  taken  by  the  sample  to  reach 
the  designated  temperature  may  be  negligible  relative  to  r.  yet  may  uffect  the  com¬ 
parison  with  the  tabulated  values  [28].  Moreover.  F{x)  may  change  in  different 
regimes  of  the  «(r)  curves,  as  mentioned  earlier. 


A  new  method 

Here  we  suggest  a  step  by-step  ttpprouch  to  determine  F(x).  It  wus  noted  by 
Hancock  and  Shurp  [28]  thut  for  many  forms  of  /’(»).  the  plot  of  log]  —  Inf  I  —  a)] 
vs.  log  t  is  almost  linear  if  a  is  restricted  to  between  0,15  and  0.5.  Using  a  com¬ 
puter  program  to  generate  artificial  values  und  their  log-ln  plots,  we  have  found 
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this  I  rue  lor  all  t  lie  theoretical  forms  listed  in  Table  I.  unit  the  exception  of  H. 
C  and  F.  The  slope  in  each  ease  is  listed  there  under  the  C  oliimn  'in'. 

Obviously  a  log-ln  plot  is  not  very  sensitive.  II'  we  were  to  rely  solely  on  it  to 
discriminate  between  the  functional  forms  of  1(2).  the  experimental  daltt  would 
have  to  be  of  the  highest  quality,  A  slightly  more  sensitive  way  is  to  plot  [  —  I n ( I  - 
-  a)]1 1,1  vs.  i.  but  then  m  can  only  be  obtained  iteratively.  Fortunately,  m  does 
differ  significantly  between  different  groups  of  Fix!,  und  the  final  discrimination 
is  easily  achieved  by  a  further  graphical  step.  There  are  three  possible  situations 
lor  this  second  step:  ~ 

1.  hi  f  2  or  log. In  concave  upwards 

It  will  be  seen  from  Table  I  that  this  situation  suggests  an  autocatulviic  reuclion 
lor  which  i  -  A'xn(l  -  x)'1  for  certain  />  and  q.  From  the  experiment!!1  cluiii  of  j! 
ami  x.  one  can  then  do  a  leaM-squarcs  lit  on  the  graph  of  I  log  *  ,1  |t)g  x  against 
J  log(l  -  x)/ A  log  a.  and  find  p  from  the  i-mtercept  and  </  from  the  slope.  Here 
J  log  i  m  log  i(/,)  -  log  i(i-i).  etc. 

For  an  i  expression  of  this  form,  one  has  /;( I  ->')*=  r/a'.  where  a'  is  the  value 
ut  maximum  1  Using  this  relation  to  reduce  the  number  or  unknown  parameters 
to  one.  one  cun  use  a  simpler  graph  to  determine  p  and  q  [26],  However,  the  cal* 
culated  values  of  p  and  q  are  then  subject  to  the  accuracy  of  a'  and.  more  fundu- 
mentully.  the  possibility  that  p  and  q  may  change  from  one  range  of  a  to  unother 
is  not  allowed  for.  As  mentioned  above,  those  types  of ,t\\  -  a)  in  which  q  ■*  0 
represent  the  acceleratory  period  which,  in  general  should  be  followed  by  a  decay 
period  governed  by  a  different  form  of /U  -  x), 

2.  m  &  1 

The  reaction  is  either  phase-boundary  controlled  or  unimoleculur.  und  i  =■ 
«*  A(l  —  a)",  as  seen  in  Table  I.  One  then  draws  the  graph  of  log  i  against 
log(l  -  a)  to  find  /•,  the  upparent  reaction  order. 

3.  in  a  0,5 

The  reaction  is  diffusion  controlled  (see  Tuble  I).  One  has  to  test  separately 
whether  x(t  I  is  parabolic  (the  diffusion  is  in  one  dimension  I.  or  -i  ln(l  -  x)  - 
=»=  Ad  -  a)*  with  .v  =  0  (two  dimensions)  or  a  |  3  (three  dimensions). 

The  correlation  coefficient  in  the  least-squares  lit  serves,  by  measuring  the  li¬ 
nearity  of  the  i  graph,  as  a  quantitative  indication  of  the  confidence  to  be  attached 
to  the  identified  form  of/d  -  x).  It  may  be  that  p  and  </,  r.  or  s  change  once  or 
twice  us  the  reaction  proceeds  from  beginning  to  completion,  but  the  i  graphs 
will  show  it  by  displaying  several  linear  segments.  If  however,  a  part  of  the  graph 
say  from  (a(.  /, )  to  (x2.  r2)  is  non-linear,  then  Ft  j)  has  changed  to  u  form  in  another 
group,  The  first  step  should  then  he  repeated  for  iltat  part :  log  [  —  Inf  I  -  [x  s-  x,  ] ' 
/(I  -  >.,])]  is  plotted  r.v.  log(r  -  ),)  for  x  between  x,  +  O.I5(Xj  -  Xj)  und  x,  +0.5 
(x.  -  x,).  followed  by  one  of  the  above  three  alternative  procedures.  On  the  other 
hand,  if  a  good  fit  is  found  with  values  of />  and  q.  r.  or  v  that  tire  no)  in  Table  1. 
the  experimenter  should  assess  whether  theoretical  justification  cun  be  provided. 
In  this  way  new  rate  laws  may  be  identified. 

Confirmation  is  carried  out  by  plotting  the  selected  functional  form  or  forms 
on  top  of  the  experimental  curve.  A  slight  misfit  in  the  very  early  part  (~  1 
minute,  depending  on  the  sample  si/e  und  the  environment!  may  he  aitributcd 
to  thermal  lag-time  und  ignored.  The  graphs  give  K.  and  from  a  set  of  A  values 
at  different  temperatures  A ,  and  E  can  he  determined.  Note  that  the  determination 
of  these  macro-kinetic  constants  depends  on  the  form  of  /)  I  -  x)  chosen,  as  it 
should  he. 
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It  may  he  added  that  we  see  it  an  immediate  possibility  to  have  full  automation 
in  the  acquisition  and  processing  of  data  in  thermul  analysis  experiments,  The 
hardware  can  be  under  the  control  of  microprocessors  or  dedicated  minicomputers, 
and  their  output  would  go  into  a  computer  or  the  same  minicomputer.  A  com¬ 
puter  program  can  then  reduce  the  data  to  x(f)  or  id)  curves,  and  further  analyse 
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Fig.  1.  Decomposition  of  a  single  crystal  of  AgN|  at  551  K;  I.  Thermogravlmetrlc  data. 
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Fig,  2.  Log-In  plot 


the  curves  to  identify /{l  -  a)  and  so  calculate  K *  and  £.  according  to  the  method 
proposed  here,  Nevertheless,  the  physical  interpretation  of  these  results  by  the 
experimenter  remains  the  crucial  step. 

The  method  has  been  applied  to  investigate  the  kinetics  of  slow  thermal  de¬ 
composition  in  silver  azide  single  crystals  using  thermogravimetric  data,  Curve  1 
in  Fig.  t  is  a  typical  experimental  curve  of  weight  loss  against  time  /,  and  Curve  2 
is  the  corresponding  reduced-time  plot  of  x  vs,  t; r.  In  Fig.  2  we  plot  log[-ln(l  - 
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Fig.  3.  i  graph 


-  a)]  vs.  log(f/r)  for  0.15  <  x  <  0.5,  obtaining  essentially  a  linear  graph  of  slope 
*1.1.  The  decomposition  of  AgNa  therefore  appears  to  be  phase-boundary 
controlled  or  unimoleculur.  Accordingly,  we  draw  in  Fig.  3  the  graph  of  log  i  vs, 
log(l  -  a),  It  shows  that  between  I  -  a  =■  I  to  0.1.  d  =  A(!  -  *)!  with  the 
correlation  coefficient  among  the  data  points  better  than  +0.9.  Indeed,  the  plot 
of  1  -  (1  -  a)1  vs.  i  (Curve  3,  Fig.  1)  is  a  good  straight  line,  with  u  correlation 
coefficient  of  +0,99  for  0.1  <  x  <  0.9,  The  decomposition  is  thus  of  the  contract¬ 
ing  cylinder  type.  We  have  also  obtained  K ,  and  E.  and  found  that  their  inter¬ 
pretation  in  terms  of  F.qs  (3)  and  (4)  leads  to  a  plausible  physical  picture.  Further 
details  are  given  in  a  separate  paper  devoted  to  the  mechunism  of  the  thermal  de¬ 
composition  of  AgNa  (30], 

In  a  second  paper  [29],  we  discuss  the  analysis  of  dynamic  kinetic  data  which 
arc  easily  obtained  from  modern  thcrmounalytical  instruments. 


Conclusion 

We  have  described  a  systematic  method  for  determining  the  kinetics  of  solid- 
state  reactions  from  isothermal  data.  It  may  be  stressed  once  more  that  whenever 
possible  a  judgement  should  then  be  made  on  the  consistency  of  the  Km  gno  £ 
values  with  the  implication  of  J\\  -  a)  regarding  the  likely  mechanism  of  the 
reaction. 
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ANALYSIS  OF  DYNAMIC  KINETIC  DATA  FROM  SOLID-STATE 

REACTIONS 
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The  kinetics  of  heterogeneous  reactions,  involving  one  reactant  in  the  solid  phase, 
usually  follow  the  law  i  m  K*  e*p(-£/fc7’)/(l  -  j),  where  a  is  the  degree  of  conversion 
of  the  solid,  and  and  E  are  the  kinetic  constants.  A  critical  examination  Is  given  of 
the  various  methods  which  are  currently  used  to  analyse  dynamic  experimental  data. 
The  limitations  of  these  methods  and  their  Insensitivity  to  the  form  of /(I  -  x)  are 
pointed  out.  An  alternative  approach  free  from  these  limitations  Is  suggested.  In  this. 
/(I  —  x)  is  determined  from  Isothermal  experiments,  and  then  the  dynamic  data  are 
accurately  analyzed  to  obtain  the  values  of  the  kinetic  constants.  A  case  study  is  given 
to  elucidate  the  applicability  of  the  approach. 

There  are  many  reactions  of  interest  in  which  one  of  the  reactants  is  in  the  solid 
phase.  These  reactions  can  be  classified  variously  os  decomposition,  dehydration, 
calcination,  dehydroxylation,  reduction,  polymeric  inversion  and  degradation, 
oxidation  etc.,  and  they  occur  in  a  wide  range  of  substances  including  ceramics, 
explosives  and  biological  materials.  The  dynamic  method  of  studying  their  kinetics 
involves  measuring  the  reaction  rates  under  conditions  of  a  continuous  temperature 
change, 

Originally  proposed  by  Skramovsky  [1  ],  the  dynamic  method  is  becoming  in¬ 
creasingly  popular,  especially  with  the  development  of  differential  thermoanalytica! 
techniques  like  DSC  and  DTA,  Since  the  initial  temperature  can  be  chosen  so  that 
the  reaction  rate  is  relatively  insignificant  to  begin  with,  it  does  not  suffer  from 
zero-time  inaccuracy -a  problem  which  exists  in  isothermal  experiments  where  the 
temperature  is  raised  rapidly  and  then  held  constant  at  a  particular  value,  A  further 
advantage  is  that,  provided  the  dynamic  datB  have  been  unambiguously  and  cor¬ 
rectly  analyzed,  any  changes  in  the  kinetic  constants  wi II  not  he  over-looked  even 
within  small  temperature  intervals,  In  contrast,  the  isothermul  method  only  pro¬ 
vides  values  averaged  over  discrete  points  in  temperature.  Also,  when  the  method  of 
analysis  used  is  such  that  the  kinetic  constants  are  calculated  from  each  dynamic 
curve  then  very  few  samples  are  required;  only  a  milligram  or  so  of  the  material  is 
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needed  for  its  thermal  characterization.  If  many  runs  are  indeed  carried  out.  differ- 
enccs  between  individual  samples  can  be  determined.  The  last  two  advantages  are 
particularly  useful  in  single  crystal  work. 

On  the  other  hand,  intrinsic  differences  should  be  carefully  distinguished  from 
the  effects  of  experimental  conditions.  In  the  first  case  literature  data  have  shown 

that  experimental  parameters  such  us  sample  mass  [S3]  and  shupe.  particle  size  in 
the  cuse  of  powder  samples,  and  umbient  atmosphere  (or  vacuum)  can  affect  signif- 
icuntly  the  calculated  values  of  kinetic  constants  (see  [2]).  If  this  happens,  then 
whenever  possible  the  empirical  results  should  be  extrapolated  to  refer  to  a  ‘sum 
dim!' set  of  experimental  conditions.  Secondly,  the  healing  rate  very  often  affects 
the  shape  of  the  dynamic  curve  obtained,  but  discussion  will  be  deferred  to  the  next 
section.  Lustly.  the  very  fact  that  temperature  is  now  a  variable,  in  addition  to  time, 
complicates  the  analysis  of  data,  IT  due  cure  is  not  taken,  either  inaccurate  or  totally 
misleading  values  are  obtained.  In  fuel,  a  survey  of  the  literature  reveals  several 
instances  of  high-quality  experimental  data  being  mix-interpreted  by  methods 
beyond  their  ranges  of  validity,  In  this  paper,  we  first  describe  the  various  methods 
and  point  out  their  limitations.  We  thin  defend  the  approach  in  which  use  is  made 
of  both  isothermal  and  dynamic  experiments,  The  analysis  of  isothermal  data  yields 
/( I  -  x)  (as  defined  below)  unambiguously  and,  knowing/d  -#),  one  can  calculate 
individual  values  of  the  kinetic  constants  from  each  set  of  dynamic  data.  We  may 
mention  thut.  historically,  isothermal  experiments  were  the  only  ones  employed 
in  the  pioneering  age  in  the  twenties  and  thirties,  when  solid-state  reactions  began 
to  be  studied  from  the  modern  point  of  view,  as  distinct  from  that  of  Langmuir. 
Nernst  and  Taminann. 


Kinetic  equation 

As  discussed  elsewhere  [3]  the  kinetics  of  a  reaction  proceeding  isothermully  can 
usually  be  described  by  the  empirical  relation: 

(*i)iioihirNi>i  Kt.A I ~a)exp( —ElkT) .  (IV 

Here  at  is  the  fraction  of  the  solid  reacted  ,  k  Boltzmann's  constant.  T  the  tempera¬ 
ture.  and /(I  -*)  and  the  constants  Kt  and  E  are  characteristic  to  the  reaction.. 
The  function /U  -a)  may  change  in  different  ranges  of  a  but  is,  for  a  given  a,  inde¬ 
pendent  of  T.  at  least  within  a  range  of  T.  Kr  and  E  should  be  the  same  for  the 
same./]  I  -  a),  If  the  rate  controlling  step  of  the  reaction  occurs  on  the  reactant-free 
surface  or  on  the  reactant-product  (solid)  interface,  then  KT,  will  contain  the  sur- 
fuce-to-volume  ratio.  In  other  words,  the  reacting  system  should  really  be  normal¬ 
ized  per  unit  area  rather  than  per  unit  size,  and  K ,  be  given  in  units  such  as  mole¬ 
cules  s" 1  m"11. 

Some  authors  have  questioned  the  general  validity  of  (1)  on  various  grounds 
[4-7],  However,  in  the  literature  (I)  is  almost  always  successfully  fitted  to  experi¬ 
mental  data.  Indeed,  this  empirical  relation  can  be  given  mechanismic  jutification 
(see  [3]). 

In  dynamic  experiments,  also,  it  is  commonly  agreed  that  (i)  may  be  adapted  to 
describe  the  reaction  rate : 


&  -  K.rA\-*)«*pl-Elknt)l  (2) 

nn  is  controllable  by  the  experimenter.  Some  temperature  programs  offer  the 
muthematicul  advantage  that  exp  [-E/kTyt  can  be  integrated  analytically  (see 
later  discussion  on  integral  methods  of  dulu  analysis).  Examples  are  the  hyperbolic 
program  where  IT-  A  -  Bi,  i.e.  t  •=  BTl  [8,  9];  a  parabolic  program  in  which 
AT1  +  BT  -  C  -  f.  i.e.  I/T  -  2,47"  +  B,  with  B  -  AEiK,  by  iteration  (10):  and 
an  exponential  program  so  that  t  =>  exp  ( -  B!Ti  with  B  a  E;k  [10].  For  the  sake 
of  experimental  convenience,  however,  the  arrangement  is  usually  that  t  -  <t> 
(.4.  C  and  <t>  represent  constants  in  a  particular  run  of  experiment). 


ii  ':*i.  ill ajiwn.sr', 
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However,  it  has  been  taken  by  some  authors  who  object  to  (2).  that 

da"|^)d'  +  l^)dr+(-;4!^  <3) 

/  t'3f  | 

where  ---  ■  (S)tioth«imiif  The  argument  is  then  that  (2)  is  seen  to  be  inade- 

1  < 1  'T.+ 

quate  even  in  the  case  of  dd  -  0,  since  the  second  term  on  the  R.H.S.  of  (3)  is  non* 
zero  but  left  out  in  (2):  see  e.g.  [I  I.  12).  It  has  further  been  proposed  [13],  by  a 
derivation  starting  from  (3),  that  (2)  is  correct  only  if  it  includes  the  extra  fuctor 

[1  +  (1  -  T,IT)ElkT) 

in  which  7„  a  7(0).  Rut  wc  hold  that  (3)  is  unsound.  Given  /.  7  und  <!>,  a  is  not 
uniquely  determined  and  therefore  not  a  function  of  these  system  variables:  dda  > 
>  0.  as  a  cannot  decrease  even  for  negative  d7  and  dd!  Nevertheless,  the  inexact 
differential  da  can  be  integrated,  if  the  dynamic  process  can  be  treated  as  the  limit* 
ing  case  of  a  series  of  time  intervals,  during  which  the  reaction  proceeds  isother* 
mally  according  to  (1)  but  at  the  end  of  each  of  which  T  is  altered,  in  a  time  so  short 
that  during  it  the  sample  is  unchanged.  Along  this  path  P  the  result  is  easily  ob¬ 
tained  [14]: 

«  r 

«  p  r. 

where  K  m  K*.  exp (~E/kT).  We  must  emphasize  that  (4)  is  not  logically  self* 
evident,  us  is  sometimes  implied  [14]  or  argued  by  mathematical  operations  based 
on  the  presumption  that  a  -  a(7,  t)  [15].  Rather,  it  comes  from  the  assumption 
that  the  reuction  under  study  Involves  no  stow  processes,  so  that  dl  depends  only 
on  the  present  values  of  a  and  T(A  is  a  function  of  state),  but  not  on  the  history  of 
the  reacting  system  (c.f.  [16]).  Only  by  this  assumption  (absence  of  memory  effects) 
can  the  dynamic  process  be  treated  as  P,  (Experimentally,  a  temperature  program 
with  temperature  jumps,  which  approximates  P,  has  been  realized  on  a  thermo- 
balance  interactively  controlled  by  computer  [17].)  Equations  (2)  and  (4)  are  of 
course  equivalent.  Their  validity  bus  also  been  shown  by  'rational*  thermodynamic 
arguments,  in  which  the  functionul  relation  i  -  g  (r.  Rt)  is  regarded  as  the  ‘con¬ 
stitutive  equation’  characterizing  the  reaction  system  [18].  In  a  new  direction. 

the  possibility  has  been  suggested  [64]  that  solid-state  reactions  may  be  studied 
by  far  from  equilibrium  thermodinamics. 

On  the  other  hand,  experience  shows  that  apparently  (4)  is  not  always  followed 
exactly.  Consider  a  reaction  being  investigated  by  a  series  of  experiments  conducted 
at  different  heating  rates  d  but  with  the  same  initial  temperature  7],.  In  (4)  we  see 
that  the  R.H.S.,  for  a  given  upper  temperature  limit  7,  is  directly  proportional  to 
I/d.  Plots  of  the  L.H.S.  vs,  T  should  therefore  all  have  the  same  shape.  It  may  hap¬ 
pen.  however,  that  increasing  departure  from  isomorphism  is  seen  when  experimen¬ 
tal  data  obtained  at  higher  d  are  so  analyzed .  The  most  probable  explanation  is  that 
the  temperature  change  is  too  fast,  causing  the  temperature  distribution  in  the 
sample  to  become  significantly  non-uniform.  In  fact,  thermal  equilibrium  is  an 
underlying  assumption  when  (4)  is  derived  above:  without  it  i  will  depend  on  the 
thermal  history  of  the  reacting  system, 

Other  factors  may  also  be  at  work.  The  reaction  rate  may  be  sensitive  to  the 
structure  of  the  reactant,  and  a  higher  d  can  enhance  the  defect  density  or  change 
the  activation  energy  of  reaction  at  a  defect  site  (19).  In  branched-chain  reactions, 
the  speed  of  the  progressive  accumulation  of  active  centres  may  vary  with  the  rate 
of  change  in  temperature  [20].  If  the  reaction  is  a  surface  process,  the  distribution 


of  reaction  'centres*  among  corners,  edges  or  faces  of  the  sample  may  change  with 
<j>  [21],  It  may  also  be  that  the  chemical  system  under  study  has  multiple  reactions 
proceeding  concurrently  in  it.  and  they  have  different  E  (22].  All  theae  variations 
in  E  may  be  accompanied  by  changes  of  K *  in  the  same  direction.  Because  of  this 
coupling,  a  linear  relation  between  E  and  log  KT  is  sometimes  observed.  Called  the 
‘compensation  effect',  this  phenomena  does  not  necessarily  mean,  as  was  suggested 
[4],  that  the  Arrhenius  expression  in  the  R.M.S.  of  (4)  is  invalid.  In  all  these  cases, 
by  varying  <j>  the  experimenter  can,  in  fact,  gain  additional  Insights  into  the  mecha¬ 
nism  of  the  reaction,  or  distinguish  between  the  competitive  reactions  in  the  react¬ 
ing  systems  (u  situation  usually,  though  not  always,  indicated  by  the  presence  of 
multiple  peaks  in  the  d  curves).  This  is  possible  if  the  method  of  data  analysis 
employed  is  such  that  K„  and  F.  are  determined  from  a  single  d  or  *  curve,  rather 
than  from  data  at  a  number  of  heating  rates.  The  method  wc  suggest  will  be  of  this 
type. 

It  may  also  happen  that  in  (2)  K*  x  T,  so  that  It  cannot  be  taken  outside  the 
integral  sign  in  (4),  Indeed,  modifications  have  been  suggested  of  some  methods  of 
data  analysis  (those  that  assume  a  reaction  order  for  the  reaction)  to  take  this  extra 
temperature  dependence  into  consideration  [23].  However,  even  when  theoretically 
required,  the  correction  may  for  practical  purposes  be  ignored,  unless  £  is  small  or 
temperatures  used  are  very  high  (d  In  K/dT  -  [E  +  kT]lkT ').  Likewise,  any  slight 
temperature  dependence  of  E  can  usually  be  neglected.  Furthermore,  irrespective 
of  this  or  the  above  complications  the  form  of/[  I  -  a)  in  (2)  and  (4)  it  not  affected. 
It  should  be  the  same  as  that  in  (I),  on  the  basis  that  the  dynamic  process  can  be 
treated  as  the  limiting  case  of  a  series  of  Isothermal  intervals,  as  already  mentioned 
above. 


Data  analysis 

Experimentally,  <2  or  a  is  obtained  by  DSC,  DTA,  TO,  DTO,  quantitative  IR 
spectroscopy  or  X-ray  diffraction,  dliatometry,  or  measurements  of  chemi-lumines- 
cence,  ultrasonic  attenuation,  dielectric  constant,  viscoelasticity,  thermal  or  electri¬ 
cal  conductivities,  or  optical  reflectivity  when  changes  in  these  characteristics  can 
be  correlated  with  a,  The  oldest  technique  is  thermomanometry,  in  which  the  p-es- 
sure  of  an  evolved  gas  is  measured  at  constant  volume,  but  its  use  has  so  far  been 
more  popular  in  isothermal  experiments.  Many  methods  of  analyzing  d  or  a  data 
have  been  proposed  to  calculate  the  kinetic  constants  E  and  Km  (for  a  critical 
review  of  the  earlier  work  see  [24])  and  sometimes  also/\l  -a).  Often  they  were 
originally  formulated  with  reference  to  one  particular  instrumentation,  but  they 
may  be  made  generally  applicable  to  all  techniques  after  quantities  measured  on 
DSC,  DTA.  TG  instruments  etc.  are  all  interpreted  in  terms  of  at  and  d.  On  the 
other  hand,  their  validity  does  depend  on  the  particular  reaction  whose  data  are 
being  analyzed.  Their  limitations  in  this  respect  form  the  subject  of  our  discussion 
below.  They  will  be  examined  in  three  groups:  peak-temperature,  integral  and 
derivative  methods,  in  this  order.  Sophisticated  instrumentation  systems  are  com¬ 
ing  into  use,  that  incorporate  computers  to  establish  baselines  or  other  null  settings, 
to  carry  out  automatic  data  acquisition,  and  to  let  the  experimenter  interactively 
analyse  the  data  (e.g.  (25 1).  Such  advances  do  not.  however,  remove  the  danger  of 
uncritical  choices  of  the  method  of  data  reduction. 
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Critical  examination  of  current  methods 
Peak-temperature  method 

Kissinger  [26]  considers  reactions  of  the  type /{I  -a)  -  (1  Differentiating 
(2)  with  respect  to  t ,  and  setting  the  resulting  expression  to  zero,  he  obtains 

im(E/k)WT>)  -  exp(-5/*rm)n(l  (5) 

in  which  m.  signifies  *pjak’  quantities,  at  the  point  of  maximum  d  where  4t  ■  0. 
He  next  assumes  that «( I 1  a  I;  therefore 

*7Vccexp(-£/*r„)  (6) 

regardless  of  /i,  which  itself  may  be  calculated  from  the  shape  of  the  *  (0  curve. 
£  and  k,  on  the  other  hund,  arc  obtained  by  performing  a  series  of  experiments  at 
different  <t>.  An  aspect,  which  we  regard  as  an  inefficiency,  of  Kissinger's  method  is 
that  only  one  point  on  the  curve  it  used  although,  in  the  case  where  multiple  peaks 
occur  signifying  that  different /[  1  -  a)  and  E  govern  different  sections  of  the  curve, 
the  method  should  still  be  applicable  to  each  peak. 

There  is,  however,  an  important  limitation.  The  a  priori  condition  that/ll  -  or)  - 
-  (I -a)’  is  actually  valid  only  in  very  special  circumstances,  namely  when  the 
rate-limiting  step  of  the  reaction  it  the  inward  movement  at  a  constant  speed  of 
the  reactant-product  Interface,  where  n  is  0.  >A  nr  2/3  for  one-,  two-,  or  three- 
dimensional  movement,  respectively,  or  when  the  reaction  is  unimolecular  so  that 
«  •  1.  Even  among  these  special  cases,  the  other  approximation  that  Kissinger 
uses  is  still  conditional,  since  a(l  -*«)"' 1  ■  l  only  for  n  -  I.  When  n  is  i/i  or  2/3, 
this  expression  varies  with  approximately  as  A«m  n(l -«)/(! -xm)J  "  i  0.2 
Aam.  where  Av.m  is  the  variation  In  itself.  In  the  Appendix  we  show  that 
changes  with  <j>  In  the  general  case,  Hence,  when  an  apparent  reaction  order  exists 
and  is  V2  or  2/3,  Kissinger's  method  can  lead  to  a  systematic  deviation  in  (6)  and  thus 
generate  a  significant  but  hidden  error  in  the  calculated  £and  K, . 

If  no  apparent  reaction  order  exists,  then  it  definitely  should  not  be  used,  other¬ 
wise  an  approximately  linear  plot  from  (6)  results  in  totally  misleading  values  of 
the  kinetic  constants,  An  example  is  in  the  decomposition  of  benzenediazonium 
chloride:  It  derives  from  DTA  data  a  value  of  £  that  Is  40%  lower  than  the  nearly 
Identical  values,  obtained  by  applying  other  methods  of  analysis  to  the  data  from 
DTA  as  well  as  other  techniques  [27],  Other  examples  are  in  the  study  of  lithium 
aluminium  hydride,  where  the  Kissinger  values  are  half  of  the  isothermal  result 
[28],  in  RDX  where  it  Is  again  40%  lower  than  all  the  values  calculated  by  other 
methods  [29],  and  in  urea  nitrate,  where  it  is  30%  lower  [30]. 

Integral  methods 

The  L.H.S.  of  (4)  is  a  function  of  *  only  and  will  be  denoted  by  F(x)i  the  R.H.S. 

can  for  practical  purposes  be  equated  with  j  K>4>  d T,  since  in  experiments  T„  will 

be  such  that  reaction  velocity  is  negligible  below  it.  i.e.  Tu  4  Elk.  In  vjew  of  these 
considerations,  many  authors  have  proposed  different  methods  of  analysing  a.(T) 
data.  j- 

The  temperature  integral  j  exp(-E/kT)dr  has  no  analytical  solution.  (In  the 
b 

unusual  case  of  a  hyperbolic,  parabolic  or  exponential  temperature  program, 
on  the  other  hand,  exp  (-ElkT)/t  is  integrate.)  The  numerical  values  of  the  inte¬ 
gral  have  been  compiled  but,  being  a  function  of  both  E  and  T.  are  not  directly 
useful  unless  an  Iterative  solution  of  (4)  by  trial-and-error  is  resc-rted  to.  Such  an 

approach  has  been  advocated  by  Zsakd  [31]  who  considers  in  particular  the  cases 

% 

of /(I  -x)  -  ( I  -  a)"  with  n  -  0,  1/3,  1/2,  2/3,  1  or  2.  when  $(*>■  log  f  — 
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has  simple  analytical  expressions,  and  by  Satava  and  Skvira,  [32]  who  generalize 
the  method  slightly  by  tabulating  the  values  of  y(a<),  0  <  *  <  1,  for  some  other 
forms  of/]  I  -  at). 

For  more  efficient  approaches,  approximations  to  the  integral  are  necessary. 
Thus,  taking  the  first  two  terms  in  an  asymptotic  (</  ■  ElkT  -*  oe)  series  of 

|  exp  ( -  r)  - Coats  and  Redfern  obtain  the  linearized  relation  [33]: 

In  (mr~ ')-  A  -  Ejk(T-')  (7) 

where  A  ■  ln(X  ,k/Ed>)(  I  -2AT/£)is ‘sensibly  constunt'  if  the  range  of  temperature 
AT  is  small.  They  further  assume  that  /]l -#)  -  (I  and  so  /]*)  can  be  cal* 
culated  at  each  (*,T).  Plotting  (7)  for  several  values  of  a  thus  gives  E  und  K,., 
making  use  of  only  one  set  of  dutn  corresponding  to  u  single  i>. 

Several  cautionary  notes  should  again  be  made  here.  The  assumption  for/]  I  -  a) 
has  already  been  discussed.  Similar  to  the  case  of  Kissinger’s  method,  results  ob- 
tained  may  be  wrong  and  misleading  If  this  functional  form  it  not  independently 
determined  beforehand.  Thus,  in  a  study  on  the  dehydroxylation  of  kaolinite  [34], 
straight  lines  over  different  ranges  of  (7"* ')  are  given  by  (7)  for  a  whole  series  of 
values  of  n,  namely.  0,  0.5,  0.667,  a  and  2.  In  particular,  plots  using  n  -  I  and 
n  •  2  are  almost  equally ‘good’. 

Secondly,  the  accuracy  of  the  asymptotic  approximation  is  rather  low.  By  com¬ 
paring  its  values  with  tabulated  values  of  the  integral  [31.  35,  36],  we  find  its  rela¬ 
tive  errors  to  be  d/// -  20%  at  u  -  5,  5%  at  u  -  10.  and  1.5%  at  u  «  20.  Thus, for 
example,  If  £  is  I  eV,  then  for  an  accuracy  of  98  %  the  highest  temperature  reached 
in  the  experimental  run  should  not  be  more  than  600  K,  a  very  low  figure  for  most 
materials  though  it  is  higher  for  larger  £.  Additionally,  expanding  A  Into  a  power 
series  shows  that  A  A/ A  a  2/rd  TIE,  so  that  at  say  2  %  Inaccuracy  the  range  of  tempera¬ 
ture,  AT,  from  which  («,  T)  points  are  selected  should  be  less  than  I0OK  (for  £  - 
I  eV).  The  total  possible  deviations  in  the  calculated  £  and  Kr  are,  to  first  approxi¬ 
mation,  the  sum  of  the  A//I  and  AA/A.  It  certainly  is  unsatisfactory  if  they  are  large' 
and  yet  nowhere  mentioned  in  the  calculation. 

Other  approximations  to  the  temperature  integral  have  been  suggested  by  van 
Krevelen  etal.  [37]  and  by  Horowitz  and  Metzger  [38],  who  made  use  of  certain 
asymptotic  expansions  in  the  vicinity  of  Tm,  the  temperature  at  peak  reaction  rate. 
Both  have  been  shown  [39]  to  be  even  lets  accurate  than  the  Coats  and  Redfern 
approach,  and  so  will  be  left  out  in  our  discussion. 

Amongst  the  integral  methods,  the  best  is  probably  the  one  due  to  Ozawa,  which 
requires  data  at  different  d  but,  in  it/] I  -  a)  remains  completely  general.  The  ap¬ 
proximation  to  the  temperature  integral  Is :  - 

r 

J  exp(-£/4ndTs>  -J  Iff  *M  (8) 

o 

so  that  from  (4) 

logdi  +  0.457(£//t)/r,  -  log  4>,  +  0.457(£/*)T3  (9) 

where  T,  and  Tt  are  taken  at  an  arbitrary  but  identical  value  of  *  in  the  two  curves 
corresponding  to  heating  rates  and  Plotting  log  <t>  vs.  \IT  for  selected  values 
of  x  should  therefore  produce  straight  lines,  the  slopes  of  which  give  £  [40]. 

Three  comments  are  appropriate  here.  By  comparing  (8)  with  tabulated  numer¬ 
ical  values,  we  see  that  it  is  7%  out  at  u  -  10  or  T  -  i  1 70K,  and  3%  and  less  only 
for  T  <  720K  (if  E  -  I  eV).  These  errors  should  be  examined  before  Ozawa's 


method  is  applied.  Secondly,  the  method  has  been  modified  [41]  to  read,  in  place 
of  (9). 

A  in  &AT,,,  '=  0.457  Eik  (10) 

in  which  A  In  </>  ■  lnd>i  -  ln^4.  etc.,  and  m  denotes,  as  before,  peak  quantities. 
This  relation  may  be  compared  with  (6)  but  In  general  it  does  not  hold  since,  as 
shown  in  the  Appendix,  xm  varies  with  4>.  Lastly,  like  Kissinger's  method,  £  cannot 
be  determined  from  data  at  a  single  </>,  and  In  some  cases  this  may  be  a  disadvantage, 
as  discussed  before. 


Derivative  methods 

The  derivative  methods  oiler  tin  udvuntugc  over  those  described  ubove  in  invok¬ 
ing  no  mathematical  approximations.  Unfortunately,  they  use  d  data  which,  with 
present  instrumentation,  tend  to  be  of  lower  quality  whether  they  arc  obtained  by 
numerically  differentiating  the  i  data  or  arc  direct  experimental  read-outs. 

The  most  straightforward,  but  us  it  stands  relatively  inefficient,  of  the  derivative 
methods  is  to  write  (2)  us :  - 

In  (d//(  I  -  *))  ■  In  K,  -  ElkT  (II) 


and  to  substitute  different  of/(l  -  a)  until  a  linear  plot  appears  [42],  Later,  we  shall 
argue,  however,  that  even  this  labor  omnia  vincit  approach  like  all  dynamic  meth¬ 
ods  in  general,  cannot  guarantee  correct  values  of  £  and  K ,  (nor  an  unambig¬ 
uous  form  of  yU  -  a)  in  this  specific  case),  although  the  labour  it  involves  may  be 
undertaken  by  the  computer. 

The  earliest  derivative  method  is  probably  that  of  Borchardt  und  Daniels,  origi¬ 
nally  formulated  for  homogeneous  reactions  in  the  liquid  phase  [43]  but  later  extend¬ 
ed  to  solid-state  reactions  [44]  for  which  It  is  now  frequently  used.  The  method 
putsyX  I  —  at)  *  (I -a)"  into  (II),  with  n  given  a  guessed  value,  and  if  a  linear  plot 
results  then  £  and  Kr  are  obtained  from  it.  Based  on  this  method.  Hauser  and 
Field  [45]  have  developed  a  computer  procedure,  in  which  plots  are  generated  for 
a  series  of  values  of  n  incremented  at  discrete  steps,  and  the  'best'  one  is  then  se¬ 
lected  to  yield  £.  K*.  and  n,  An  attraction  of  this  method  is  that  «  can  be  readily 
selected  by  eye.  Alternatively,  since  in  this  case 


A  In  d/dln(l  -x)  »  -(Eik[AT~  1  A  In ( 1  -a)])  +  n  (12) 


a  plot  of  the  L.H.S.  vs.  the  quantity  in  the  square  brackets  at  once  gives  £  from  the 
slope  and  n  as  the  .v-intercept  [47  ],  If  constant  A  Ini.  A  Inf  I  -x).  or  AT" 1  is  selected. 
Eq.  (12)  can  be  further  simplified  [61  ].  We  have  emphasized  previously  the  falli¬ 
bility  in  presuming  such  a  convenient  form  of  /f  l -*):  Ozawa  [47]  has  commented 
on  the  possibility  that  this  procedure,  and  the  integral  method  of  Coats  and  Red- 
fern.  may  give  false  values  of  £  and  K.r  ■  in  addition,  since  (12)  Involves  the  ratios 
of  differences,  the  quality  of  data  called  for  is  even  higher  than  that  demanded 
alone  by  the  use  of  i:  experimental  data  so  plotted  more  often  than  not  show  very 
large  scatter.  The  Rogers  and  Morris  method  [48]  plots  din*  vs.  T*'.and  can  be 


seen  to  be  the  speciul  case  of  n  «■  0  in  (12).  An  example  of  the  general  danger  that 
very  linear  plots  may  sometimes  appeur  even  If  the  applied  method  is  not  valid  is 
given  by  Patel  and  Chaudhri.  The  Rogers  and  Morris  method  was  used  to  analyze 
DSC  duta  on  lead  azide,  and  a  straight  line  results  although  the  calculated  £  turns 
out  to  be  1 80 u „  larger  than  the  Ozawa  value  [49).  Conversely,  the  coincidence  of 
values  calculated  by  various  methods  need  not  prove  that  these  methods  are  all 
applicable  to  the  case  in  hand.  A  counter-example  is  provided  by  a  DSC  study  on 
RDX  [29],  where  the  Rogers  und  Morris  value  agrees  well  with  other  values  but 
the  complex  decomposition  is  beyond  doubt  far  from  the  n  -  0  type. 


On  the  other  hand.  D4vid  and  Zelenyiinszki  [50]  plot  In 


d 

dr 


(I -*)  /(! 


-1 


against  (f  ');  this  amounts  to  assuming  a  'reaction  order'  »  ■  I.  It  serves  as  yet 
another  example  of  the  futility  of  linear  plots,  for  their  method  gives  such  plots  for 
the  decomposition  of  ’a  wide  range  of  materials’  including  calcium  oxalate  and 
polyethylene  which,  most  likely,  are  not  of  first  or  any  other  'order-. 


Some  of  the  integral  and  derivative  methods  described  in  the  foregoing  have  been 
compared  by  testing  their  accuracies  on  synthesised  DTA  data  (exact  as  well  as 
with  artificial  random  error)  for  one  E  value  and  temperature  range,  the  reaction 
considered  being  of  the  type  with  a  reaction  order  [51).  Among  the  methods  not 
included  there  is  that  due  to  Friedmann  [52],  It  probably  is  the  most  general  among 
the  derivative  methods.  Like  Ozawa's  procedure,  it  makes  no  assumption  about 
/( 1  —  at),  although  it  requires  d  data  which,  furthermore,  have  to  be  at  a  number  of 
Once  again,  from  (2)  with  dx/d T  ■ 

ln(ard>)  -  ln(A„/(l-a))  -  ElkT.  (13) 

Since  K,.f(  \  -*)  is  identical  for  the  same  value  of  a.  taking  ar  and  the  correspond¬ 
ing  7"  from  several  <f>  one  can  determiner. 

It  is  our  contention  that  even  Friedmann's  method  hus one  basic  limitation  which, 
more  significantly,  is  shared  by  all  dynamic  methods  described  above.  The  point 
in  question  is  that  all  of  them  have  to  presume  the  constancy  of/"(  I  -*)as  the  tem¬ 
perature  is  changed.  However,  since  mechunisms  of  solid-state  reactions  arc  gener¬ 
ally  complicated,  there  is  no  gcnerul  justification  for  this  presumption,  though  It 
may  be  true  for  particular  reactions  within  specific  temperature  range*.  An  illustra¬ 
tion  is  the  case  where  parallel  reaction  paths  exist,  each  with  values  of  AT*  and  E 
such  that  a  quantitative  change  of  T  will  lead  to  a  qualitative  change  in  the  dominat¬ 
ing  path.  Another  case  is  where  the  identification  of  the  rate-limiting  step  depends 
on  T.  Methods  have  been  proposed  which,  by  the  use  of  computers,  try  different 
forms  of /X 1  -#)  in  analyzing  the  dynamic  data  [53-55],  However,  the  search  is 
limited  to  functional  forms  which  are  already  known, 

More  importantly,  from  our  own  experience  with  azides  we  have  strong  doubts 
as  to  the  exactness  in  determinlng/f  l  -  at)  or  even  its  constancy  from  dynamic  data. 
Likewise,  in  a  study  on  the  dehydration  of  manganese  formate  [54]  for  instance, 
no  unique  form  of/Tl  -a)  and  correspondingly  no  unique  values  of  £  are  identified 
even  over  appropriately  restricted  ranges  of #,  the  criterion  used  being  minimum 
standard  deviation  in  the  Arrhenius  plot.  Further  examples  are  the  thermal  dehy- 
droxylations  of  kaolinite  [34]  and  of  magnesium  hydroxide  [56].  We  suggest  that, 
in  dynamic  experiments  since  data  are  collected  under  variable  temperature  condi¬ 
tions,  the  change  due  to  /(I  -  *)  is  inherently  masked  by  that  due  to  £(7').Thls  pitfhll 
is  illustrated  In  the  Figures.  Figure  1  a  shows  the  graphs  of  #(f)  and  its  derivative 
which  are  generated  artiflcally  according  to  the  theoretical  equation  <xT  -  AT. 


Fig,  I*.  Artificial  data  *  —  Hi.  and  corresponding  «r  ■  d»/drdata,  plotted  against  T  which 

rises  linearly  with  t 

Fig.  lb.  Arrhenius  plots  of:  (1)  *r/J(l  —  *)*'».  (2)  *r/2(l  -  e)1'1,  C3)  jr,  and  [(4)  *r/ 3a */* 
for  the  data  shown  In  Fig  Is 
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exp  i  -  E/kT)  throughout  from  *  ■>  0  to  *  «■  I .  The  values  chosen  for  K„  and  E are 
10“  und  leVrespectively.  Let  us  now  examine  how  the  data  generated  according  to 
this  relution,  which  is  of  the  type  %  —  Kt.  will  be  fitted  by  different  kinetic  equations, 
one  of  them  being  the  correct  one.  In  Fig.  lb  we  plot  against  1  G'/T  the  natural 
logarithms  of  the  following  expressions :  - 

(1)  xr/3(l -x)a'\  i.e.  assuming  I  -  (I ■  Kr.  reaction  controlled  by 
three-dimensional  contraction  of  phase  boundary: 

(2)  atT/2( I i.e.  I  -  (1 -*)'■■  -  K/type: 

(3)  xr:  the  original  assumption:  and 

14)  x'flx*’1'.  i.c.  a  ■  (Ktfi  reuction  controlled  by  e.g.  three-dimensional 

growth  of  existing  nuclei. 

Ii  is  seen  that  the  incurred  /( I  -*)  in  (I)  und  (2)  still  give  virtually  linear  plots,  with 
slightly  dilferenl  slopes:  interestingly,  curve  (4)  is  so  misleading  as  to  show  two 
"linear"  segments  with  a  seemingly  significant  transition  in  between.  Experimen¬ 
tally,  Guarini  e/a/.  (57] have  noted  that  it  is  impossible  to  ascertain  from  their  DSC 
data  whether  the  monomerization  of  9-Me-IO-AcAD  hat  an  apparent  reaction 
order  of  1, 0.67,  or  0.5,  in  all  of  which  cases  E  has  about  the  same  derived  value. 

The  suggested  approach 

In  the  foregoing  sections,  we  haveditcutsed  the  limitations  regarding  the  applica¬ 
bility  of  various  methods  that  have  been  used  to  analyse  dynamic  data.  In  many 
published  works  we  find  that  often  a  number  of  apparently  different  methods  are 
used  to  analyse  the  same  set  of  data.  However,  we  think  that  in  many  cases  this 
procedure  is  of  no  real  significance,  when  some  of  the  methods  used  are  mathemat¬ 
ically  equivalent  and  therefore  lead  to  the  same  results,  or  when  some  are  invalid 
in  the  given  situation  and  thus  lead  to  doubtful  values.  The  limitations  of  the  meth¬ 
ods  express  themselves  both  as  discrepancies  in  the  calculated  values  of  the  kinetic 
constants,  and  sometimes  as  fortuitous  agreements  when  some  of  the  methods  are 
certainly  inapplicable.  (An  extreme  example  of  the  second  situation  is  that,  for 
RDX.  the  Kissinger  value  [29]  of  E  is  near  to  that  obtained  [58]  by  plotting 
x  vs  T~ a  procedure  which  has  absolutely  no  theoretical  justification.)  Accord¬ 
ingly.  we  suggest  that  the  interpretation  of  dynamic  data  should  as  far  as  possible 
be  based  on  results  from  isothermal  experiments.  A  similar  approach  has  been 
used  for  studying  the  dehydroxylation  of  kaotinile  by  Achar.  Brindley  and  Sharp 
[63]. 

One  can  unambiguously  determine ./X I  -x)  over  the  whole  range  of  x  and  over 
the  relevant  temperature  range,  from  the  independent  analysis  of  Individual  iso¬ 
thermal  curves.  A  systematic  method  of  efficiently  Implementing  this  identification 
hus  been  proposed  by  us  [3].  It  may  also  be  noted  that  thermoanalytical  equipments 
are  equally  applicable  In  isothermal  experiments  (see  e  g.  [62])  though  they  are 
more  often  used  in  the  dynamic  mode.  The  identified  formO)  of /(I  -a)  can  then 
be  substituted  into  either  (2)  or  (4).  In  this  way,  form  the  dynamic  x  or  d  data  one 
can  then  determine  accurately  the  non-average  and  single-sample  values  of  E  and 
K, :  advantages  which  have  been  mentioned  in  the  introduction  to  this  paper. 
Moreover,  the  values  will  correspond  individually  to  different  heating  rates. 

We  applied  this  approach  to  the  spinel  formation  ZnO  +  CraO,  -*  ZnCr,04. 
A  DTA  curve  (experimental  atmosphere:  N,  at  300  mm  mercury)  was  published 
in  ishii  et  at.  [59],  who  have  also  monitored  a(/)  by  chemical  analysis  when  the 
reaction  proceeded  Isothermally  in  nitrogen  flowing  at  50  ml/min,  and  showed 
that  the  Isothermal  data  fit  [I  -  <1  -  a)"*]*  ■  Kt.  We  have  measured  K  from  the 


Table  I 


IsothcrniBl  data 

r. 

*■ 

In  K 

Uag.  C 

ritativf  units 

BOO 

1 

0 

900 

4 

1.4 

1000 

2.1 

3.1 

Table  2 
Dynamic  data 


i. 

i 

t 

lii  |  th  1 

4...  t 

1  1 

i 

ftilullvc  unlit. 

A  II 

700 

0.1 

2 

-2.6 

820 

0.13 

4 

-1.4 

900 

0.2 

7 

-0,34 

1000 

0.13 

12 

+0.76 

experimental  data  points  at  t  -  20  min  in  the  published  isothermal  plots.  From 
these  values  of  K .  given  in  Table  1,  we  calculate  a  value  of  1.5  eV  for  £. 

In  Table  ?  the  values  d  were  measured  from  the  published  DTA  curve  whose 
heating  rate  was  unspecified,  and  the  x  values  were  read  off  from  the  (a,  T)  graph 
which  Ishii  eial.  have  drawn  presumably  by  Integration.  Now,  from  their  analysis 
of  the  isothermal  data  the  governing  kinetic  equation  Is,  in  differential  form,  oi  - 
■  K/[(l  — x)-'18  -  (1 -«)-•'»].  at  least  within  the  ranges  800- 1000°  and  x  from 
0  to  at  0.6  corresponding  to  Ki  •>  0  to  a  0.07,  The  Arrhenius  plot  of  d  (I  -x)* 1,11 
I(l-x)*,;!l  -  I],  for  the  four  data  points  shown  in  Table  2,  is  indeed  a  good 
straight  line.  From  the  plot  we  obtain  £  -  1.3  eV.  In  view  of  the  probable  experi¬ 
mental  errors  and  inaccuracies  in  obtaining  data  from  the  published  graphs,  we  con¬ 
sider  satisfactory  the  reasonable  agreement  between  this  value  and  the  one  calculat¬ 
ed  from  the  isothermal  data. 


Conclusion 

Most  of  the  commonly  used  methods  of  analyzing  dynamic  data  have  been  shown 
to  be  applicable  only  under  particular  conditions.  It  has  been  pointed  out  that  to 
use  these  methods  without  considerations  of  the  range  of  their  validity  can  give 
misleading  values  of  the  kinetic  parameters.  An  approach  has  been  advocated  in 
which  use  is  made  of  both  the  dynamic  and  isothermal  data:  the  functional  form 
./HI  -x)  is  determined  from  the  isothermal  experiments.  This  form  in  conjunction 
with  the  dynamic  data,  gives  the  values  of  the  kinetic  constants. 
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Appendix 


We  wish  to  predict  how  the  value  of  a  at  peak  reaction  rate  varies  with  the  heat¬ 
ing  rate  0.  At  a  =  from  (4) 

*m  T  m 

1  -  A;  I  dr  exp  (-f/A-n  h  f<*„.  «p.  rj  -  o 

I  -  x)  0  J 


.>  /(»■ 
0 


and  from  the  fact  that  i  -  0 


rvi  ■»  i  E  $  G*p(ElkTm)  ,  _  . 

J  (I  .  "yV  ,,  ■  G(ctm,  0*  » *,)  “  0 . 

Solving  the  simultaneous  equations  df  -  0  and  dG  »  0.  we  find 


tic?,,, 

I  «V  r7G 

.V 

r(  i  , 

f  *  F  Mi 

,y- 

fr;  | 

<Stl> 

1  <'<!>  M'm 

'  fTm 

-  A  1 

c0  ; 

'  ftr„ 

arm  1 

(I4> 

(15) 

(16) 


Defining  dimensionless  quantities  V  ■  E!kTm.  fim0IK,Tm.  and  H<xm)  i 
da//l  I  -  *),  we  have  the  following : 
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r'F 
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.A'l-xJ  ’ 
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Kr 
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Hence  1*"*  V-  *(2  +  I 

d0  0  A  (7(2  +  t/).iA  1 1  aj  -  /*( 1  -  x„,)/A 


exp((7),  (17) 

(l»> 


where  A  ■  u  exp  ( V ).  Incidentally,  dTJd<t>  cr.n  to  derived  in  a  similar  wav 
The  only  case  we  find  reported  in  the  literature,  in  which  otm  is  apparen,.  Inde¬ 
pendent  of  0.  is  the  primary  recrystallization  cf  pre-ccrr, pressed  copper  (60  j.  where 
xm  ■  0.5.  In  all  other  cases,  experiments  give  changing  vm.  We  have  made  a  rough 
check  on  (18)  by  taking  the  case  of  the  decomposition  of  the  explosive  RDX  [29], 
for  which  the  Rogers  and  Morris  method  aives  E  -  2. 10  eV  and  K,  *  lO,Mi“‘. 
The  reaction  is  complex,  but  these  representative  values  are  chosen  because  they 
correspond  to  an  assumed  kinetic  equation  in  which ,f\\ -a)  -  I.  We  thus  have 
very  simply  /(am)  »  a„and/'(l  -am)  «  0.  For  0  *•  0.167  K  s'1.  amis  given  as  0.62 
and  T,„  us  5I2K:  our  calculation  shows  dxm/d0  a  -0.1  s  K'\  a  value  which 
compares  well  with  the  experimental  indication  that  dxm/d0  *  (0.60  -  0.62)/ 
,(0.333  -  0.167] sK-1. 
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